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SUMMARY 

The  overall  objectives  of  this  work  are  to  (1)  refine  and  update  an  existing  soot  formation 
model  and  (2)  incorporate  this  soot  model  into  a  code  describing  a  laminar,  opposed  jet 
diffusion  flame.  To  help  advance  the  understanding  of  chemical  limitations  of  PAH 
formation  in  the  soot  model,  phenylacetylene  has  been  pyrolyzed  in  a  single-pulse  shock 
tube.  The  formation  of  PAH's  from  aromatic-aromatic  reactions  has  been  observed  and 
evidence  has  been  obtained  for  isomerization  (equilibration)  among  polyaromatic  species. 
Thermodynamics  of  high  molecular  weight  PAH's  have  been  calculated.  Benzene 
production  and  radiation  phenomena  have  been  modeled  in  an  opposed  jet  diffusion  flame 
and  where  possible  these  calculations  have  been  compared  to  experimental  results.  A 
sectional  aerosol  model  has  been  added  to  the  diffusion  flame  code  as  well  as  related 
subroutines  on  soot  formation,  growth,  and  oxidation  and  the  debugging  of  the  code  is 
nearly  complete.  Comparisons  of  the  sectional  model  and  a  global  soot  formation  model 
have  been  made  for  premixed  flames. 

I.  INTRODUCTION 

The  overall  goal  of  this  three  year  effort  is  to  incorporate  a  soot  model  into  a  code  for  a 
laminar,  opposed  jet,  diffusion  flame  and  make  comparisons  with  experimental  data. 
Specific  objectives  in  order  to  reach  this  goal  are:  (1)  identify  reduced  kinetic  mechanisms 
which  adequately  describe  ring  formation  and  growth  processes,  (2)  obtain  and  interpret 
new  data  on  polyaromatic  hydrocarbons  (PAH)  formed  during  hydrocarbon  pyrolysis,  (3) 
refine  and  update  an  existing  soot  formation  model,  (4)  incorporate  a  radiation  model  into  a 
code  describing  a  laminar,  opposed  jet  diffusion  flame,  and  (5)  incorporate  the  soot  model 
into  the  flame  code.  In  this  second  annual  report,  progress  related  to  these  items  is  briefly 
described. 


II.  RESULTS 


II.  A.  Pyrolysis  of  Phenylacetylene 

Phenylacetylene  (0.25%  in  a  bath  of  argon)  has  been  pyrolyzed  in  a  single-pulse  shock 
tube  as  part  of  an  effort  to  understand  growth  reactions,  i.e.,  PAH  formation,  during  the 
pyrolysis  of  hydrocarbons.  Products  have  been  analyzed  using  a  gas  chromatograph 
coupled  to  a  mass  selective  detector.  Light  hydrocarbons  are  principally  acetylene, 
diacetylene,  and  benzene. 

Several  principal  product  species  have  a  mass  178  amu  and  the  most  dominant  of  these  has 
not  yet  been  identified.  The  less  dominant  species  with  this  mass  are  isomers  of 
ethynylbiphenyl,  phenanthrene,  and  methylenefluorene.  Small  concentrations  of 
napthalene,  diethenylbenzene,  naphthylacetylene,  and  acenaphthylene  are  also  observed. 
Surprisingly,  very  low  levels  of  anthracene  (a  three  ringed  aromatic  of  mass  178  amu)  are 
observed  at  low  temperatures.  This  result  supports  arguments  for  a  facile  route  to 
formation  of  phenanthrene  (a  slightly  more  stable,  three  ringed  aromatic)  in  the  presence  of 
phenylacetylene.  The  lack  of  substantial  concentrations  of  pyrene  and  larger  species 


suggests  very  rapid  reactions  between  these  species  and  phenylacetylene,  causing 
additional  growth. 


II.  B.  Formation  Mechanisms  of  Polyaromatic  Hydrocarbons 

The  preferential  formation  of  phenanthrene  during  the  pyrolysis  of  phenylacetylene  was 
unusual.  Alternatively,  during  the  pyrolysis  of  toluene,  anthracene  production  appears  to 
be  favored  at  low  temperatures,  yet  fairly  constant  ratios  of  phenanthrene  to  anthracene 
were  observed  at  elevated  temperatures.  In  Fig.  1,  the  phenanthrene/ anthracene  ratio  is 
plotted  as  a  function  of  initial  post-shock  temperature  for  both  the  toluene  and  the 
phenylacetylene  systems.  For  reference,  ratios  of  the  equilibrium  concentrations  of  these 
species  is  also  provided.  This  curve  clearly  shows  the  difference  in  the  ability  of  the  two 
single-ringed  aromatic  species  to  form  three-ringed  polyaromatics.  Mechanisms  for  the 
facile  formation  of  the  polyaromatics  are  suggested  in  Fig.  2.  In  the  case  of 
phenylacetylene  as  a  parent,  it  is  assumed  that  the  phenyl  radical  is  formed  from  H-atom 
addition  to  the  alkyl  group,  followed  by  acetylene  elimination.  Alternatively,  phenanthrene 
could  be  formed  by  a  Diels- Alder  addition  of  phenylacetylene  to  benzene,  followed  by 
direct  elimination  of  molecular  hydrogen. 

The  relaxation  of  the  phenanthrene  to  anthracene  ratio  to  near  equilibrium  conditions  at 
elevated  temperatures  strongly  suggests  the  existence  of  rapid  isomerization  reactions. 
Included  in  Fig.  1  are  calculations  for  the  toluene  case,  where  anthracene  is  assumed  to  be 
initially  formed  at  a  rate  which  is  42%  faster  than  the  phenanthrene  production.  The 
isomerization 


anthracene  -»  phenanthrene 

is  assumed  to  occur  with  a  rate  constant  of  2  x  10^  exp  (-75000  cal/mole/RT)  sec-*- 
Using  the  reverse  of  this  rate,  the  relaxation  of  phenanthrene  to  anthracene  was  also 

calculated  for  the  phenylacetylene  system.  These  results  are  also  provided  in  the  figure. 

These  different  mechanisms  (formation  and  isomerization  steps)  support  a  new  hypothesis 
for  PAH  formation  under  pyrolytic  conditions  and  in  diffusion  flames.  It  is  already 
recognized  that  there  is  a  difference  between  PAH/soot  forming  processes  in  a  diffusion 
flame  and  in  the  post-flame  zone  of  a  premixed  flame.  Specifically,  mechanisms  are 
believed  to  be  fuel-type  independent  in  a  premixed  flame,  but  are  fuel-dependent  in 
diffusion  flames.  These  new  proposed  mechanisms  shed  further  light  on  this 
understanding.  For  diffusion  flames,  fuel  type  can  have  a  major  influence  on  ring 
formation  mechanisms  at  early  times  and  at  temperatures  below  about  1600K.  But  for 
longer  times  and  particularly  at  elevated  temperatures,  pyrolytic  reactions,  perhaps 
enhanced  by  oxidative  pyrolysis,  quickly  convert  the  parent  fuel  into  the  more  stable 
aromatic  species,  isomerizations  occur  and  further  ring  growth  is  governed  by  similar 
mechanisms,  independent  of  fuel  type. 

II.  C.  Thermodynamics  of  Heavy  Polyaromatic  Hydrocarbons 

As  part  of  the  effort  of  incorporating  a  soot  formation  model  into  a  flame  code,  the  effects 
of  changes  in  the  thermodynamics  of  the  average  PAH  moiety  as  PAH's  increase  in  size  is 
being  evaluated.  In  addition,  requirements  for  accurate  estimates  of  the  thermodynamics  of 
these  species  is  being  determined.  As  an  example  of  the  types  of  calculations  performed  to 
date,  heats  of  formation,  entropies,  and  heat  capacities  per  carbon  atom  are  plotted  as  a 
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Fig.  1  Phenanthrene/Anthracene  Ratio 

Pyrolysis  Data  and  Equilibrium 
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Fig.  2 


Routes  for  PAH  Formation 


2a.  Formation  of  anthracene  from  toluene 


(anthracene) 


2b.  Formation  of  phenanthrene  from  phenyiacetylene 


function  of  numbers  of  carbon  atoms  (C.A.)  in  Figs.  3-5.  Species  selected  are  composed 
of  six  ring  aromatic  structures,  and  the  heaviest  species  are  fully  symmetric,  six-sided 
polyaromatics.  Calculations  were  performed  for  these  planar  species  for  PAH's  containing 
as  many  as  864  carbon  atoms  and  72  hydrogen  atoms.  This  species  is  composed  of  397 
aromatic  rings  and  has  a  mass  of  10,440  amu.  The  asymptotic  values  at  300K  obtained 
from  Figs.  3-5  can  be  compared  with  known  values  for  graphitic  carbon.  The  heats  of 
formation  of  the  PAH's  appear  to  asymptote  near  a  value  of  1.7  kcal/mole  per  carbon  atom, 
which  is  slightly  greater  than  the  heat  of  sublimation  of  a  carbon  atom,  1.45  kcal/mole/K. 
The  limiting  entropy  is  about  2  entropy  units  (e.u.),  which  is  higher  than  the  entropy  of 
solid  carbon,  i.e.,  1.4  e.u.  Similarly,  the  limiting  heat  capacity  is  about  2.2  cal/mole/K  per 
carbon  atom,  slightly  higher  than  2.05,  the  heat  capacity  of  pure  carbon.  The  differences 
between  the  calculated  values  and  the  values  of  graphitic  carbon  are  presently  being 
evaluated. 

Typical  soot  particles  have  H /C  molar  ratios  of  about  1/8.  In  Fig.  6  are  plotted  H/C  ratios 
of  the  polyaromatics  examined  above.  Consequently,  one  might  speculate  that  soot 
thermodynamics  could  be  adequately  represented  by  the  thermodynamics  of  a  species  with 
the  appropriate  H/C  ratio,  i.e.,  C384H4g. 

Differing  thermodynamics  can  have  a  substantial  effect  on  the  equilibrium  situation  between 
a  fuel,  molecular  hydrogen  and  the  polyaromatic  product.  For  example,  consider  a  30% 
methane  mixture  in  70%  nitrogen,  initially  at  1500K  and  1  atmosphere.  Equilibrium 
temperatures  and  hydrogen  concentrations  are  plotted  in  Figs.  7  and  8  as  a  function  of  the 
polyaromatic  assumed  to  be  in  equilibrium  with  the  methane  and  hydrogen.  As  can  be  seen 
from  these  figures,  temperatures  can  vary  more  than  150K  and  hydrogen  concentrations 
can  differ  by  nearly  a  factor  of  two  for  different  assumed  polyaromatics. 

n.  D.  Benzene  Predictions  in  Methane/ Air  Diffusion  Flames 

Last  year  it  was  shown  that  benzene  formation  during  the  oxidative  pyrolysis  of  methane 
was  due  to  reactions  between  C3HX  species.  Using  this  result  a  simplified  benzene 

production  model  was  constructed  *  together  with  several  steps  describing  benzene 
pyrolysis  and  oxidation.  This  mechanism  was  added  to  a  methane  kinetics  set  (including 
C2*hydrocarbon  chemistry)  which  has  been  used  previously  for  modeling  methane 
opposed  jet  flames.  As  described  in  Appendix  A,  benzene  profiles  were  then  predicted  in 
an  opposed  jet  diffusion  flame.  Peak  benzene  profiles  are  shown  in  Fig.  9  for  flames  with 
very  low  stretch  through  extinction.  These  profiles  are  plotted  as  a  function  of  J,  which  is 
a  parameter  nearly  proportional  to  strain  rate.  As  shown  in  Fig.  9,  as  flame  stretch 
increases,  benzene  production  initially  drops  rapidly  and  then  decreases  more  slowly  with 
increasing  stretch.  Du,  Axelbaum  and  Law^  found  that  the  sooting  limit  occurs  in  the 
region  of  rapid  benzene  decrease  with  flame  stretch,  and  these  results  are  therefore 
qualitatively  consistent  with  linkage  between  benzene  production  and  soot  limits.  The 
sooting  limits  as  reported  by  Du,  et.  al,  are  reproduced  in  Fig.  9. 

Temperatures  at  the  location  of  peak  benzene  concentrations  are  plotted  along  with  peak 
flame  temperatures  in  Fig.  10  for  flames  of  varying  stretch.  This  figure  indicates  that 
benzene  production  in  diffusion  flames  occurs  at  temperatures  near  1400K,  much  lower 
than  1800K  at  which  benzene  was  formed  in  the  shock  tube  experiments  on  rich  methane 
oxidation.  Since  dominant  kinetic  pathways  can  be  altered  by  a  substantial  shift  in 
temperature,  the  full  kinetics  set  (all  known  mechanisms  of  benzene  formation  included) 
were  used  to  recalculate  selected  diffusion  flames.  Under  all  conditions  only  trivial 
changes  were  observed  in  the  benzene  profile,  confirming  the  earlier  result  of  the 
importance  of  propargyl  species  to  benzene  production  in  methane  systems. 
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Fig.  3  AHf/C.A.  vs  Carbon  Atoms 


Fig.  4  Entropy/C. A.  vs.  Carbon  Atoms 
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.  Fig.  5  Heat  Capacity/CA  vs.  Carbon  Atoms 
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Rg.  6  H/C  Molar  Ratio  vs.  Carbon  Atoms 
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Fig.  7  Dependency  Of  Equilibrium  Temperatures 
On  PAH  Thermodynamics 


Fig.  8  Dependency  of  Equilibrium  Hydrogen 
On  PAH  Thermodynamics 
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Fig.  9  Peak  Benzene  As  Function  Of  Strain 
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Fig.  10  Temperature  vs.  Strain  Rate 
Laminar  Opposed  Jet  Diffusion  Flame 
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In  addition,  an  attempt  was  made  to  model  a  nitrogen  diluted  methane  studied  by  Zhang, 
Atreya  and  Lee^.  All  major  species  were  predicted  well.  As  shown  in  Fig.  1 1,  shapes  of 
C2-species  were  described  well,  but  typically  peak  concentrations  were  modeled  only 
within  a  factor  of  two  to  three,  perhaps  due  to  the  use  of  only  a  skeleton  C2-hydrocarbon 
mechanism.  Benzene  profiles  were  predicted  very  well  (shape  and  peak  concentration) 
although  this  result  may  be  somewhat  fortuitous  due  to  the  inadequacy  in  the  comparison  of 
the  C2-species. 

In  general,  these  kinetic  results  are  very  encouraging  for  the  soot  modeling  effort  in  that 
opposed  jet  diffusion  flames  can  be  modeled  fairly  inexpensively  with  chemistry  describing 
species  as  large  as  benzene.  For  example,  in  order  to  determine  benzene  concentrations, 
etc.,  as  a  function  of  strain  rate,  567  flames  were  calculated  in  12  hours  using  an  IBM 
RS6000  workstation. 

II.  E.  Radiative  Transfer  in  Sooting  Opposed  Jet  Flames 

Our  previous  investigations  of  radiative  transfer  in  non-sooting  opposed  jet  flames  showed 
that,  for  realistic  values  of  strain  rate,  up  to  5%  of  the  flame  enthalpy  release  could  be 
converted  to  thermal  radiation.  The  accompanying  reduction  in  peak  flame  temperatures 
gave  rise  to  substantial  reductions  of  about  33%  in  the  peak  NO  concentrations^-  At  very 
low  strain  rates  and  high  pressures,  where  the  pressure-flame  width  product  can  be  one 
atm -cm  or  more,  it  was  shown  that  optical  thickness  effects  can  become  important  and  an 
expression  was  derived  for  self-absorption^.  Radiative  loss  effects  can  be  particularly  large 
in  sooting  flames.  To  assess  the  effect  in  sooting  opposed  jet  flames,  the  gas  band 
expressions  for  the  net  cooling  rate  have  been  extended  to  include  soot.  The  problem  has 
traditionally  been  a  difficult  one  because  of  the  overlaps  between  the  broadband  soot 
absorption  profile  and  the  molecular  resonances.  Net  gas  cooling  results  from  radiative 
emission  minus  self-absoprtion.  The  self-absorption  term  has  been  derived  from  a  solution 
for  the  radiative  intensity  by  taking  hemispherical  and  wavelength  averages.  The  result  is  a 
semi-analytical  expression  for  the  net  cooling  rate  in  which  the  emission  is  offset  by  soot 
and  gas  self-absorption,  and  by  soot-gas  radiative  interchan ge®.  Figure  12  shows  the 
calculated  radiative  loss  in  a  model  sooting  diffusion  flame  in  which  a  band  of  soot  has 
been  synthetically  inserted  between  the  flame  and  the  stagnation  plane,  and  the  peak  soot 
volume  fraction  varied  parametrically.  The  flame  is  at  10.5  atmospheres  and  the  strain  rate 
is  20  sec'*;  the  fuel  and  oxidizer  temperatures  have  been  chosen  to  give  a  peak  flame 
temperature  approximating  that  of  jet  fuel^.  The  assumed  location  of  the  sooting  region  is 
such  as  to  produce  an  effective  radiation  temperature  for  the  soot  of  about  1500K,  about 
900K  cooler  than  the  peak  flame  temperature.  Gas  band  self- absorption  is  important  at  this 
combination  of  pressure  and  strain  rate.  For  peak  soot  volume  fractions  in  excess  of 
roughly  lO-^  SOot  radiative  transfer  begins  to  dominate,  and  predicted  radiative  loss 
fractions  rise  to  the  level  of  20%.  At  the  highest  volume  fractions,  the  reduction  in 
temperature  in  the  sooting  region  is  large  enough  to  raise  the  possibility  of  strong 
interactions  between  soot  growth  and  soot  radiation.  These  effects  are,  hov'ever,  strongly 
sensitive  to  the  presumed  location  and  radiative  temperature  of  the  soot.  A  self-consistent 
analysis  will  be  provided  by  the  coupled  soot  growth-opposed  jet  flow  code  whose 
development  is  described  in  Section  II.G. 

II.  F.  Comparison  of  Soot  Growth  Models 

The  soot  growth  model  developed  under  this  program  combines  a  simplified  picture  of 
particle  inception  linked  to  local  benzene  production  rates  and  an  aerosol  dynamics-based 
treatment  of  particle  growth  using  the  algorithms  of  the  well-known  MAEROS  code 
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Fig.  11  Comparison  of  Model  and  Experiment 
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This  approach  necessitates  solving  a  dynamical  equation  for  each  particle  size  class  or 
section.  Other  investigators  have  proposed  simpler,  global  models  based  on  assuming  a 
monodisperse  size  distribution  with  two  dynamical  equations,  one  for  the  particle  number 
density,  and  the  other  for  the  volume  fraction.  One  such  model  is  that  due  to  Linds  tedA 
which  features  an  inception  rate  proportional  to  the  acetylene  concentration,  and  a  surface 
growth  rate  proportional  to  the  square  root  of  surface  area.  It  has  given  good  agreement 
with  sooting  opposed  jet  data,  and  has  been  incorporated  into  laminar  flamelet/p.d.f. 
models  for  turbulent  jets^.  Test  comparisons  of  this  model  and  the  sectional  model  have 
been  made  for  the  premixed  flame  data^-lO  The  first  comparison,  shown  in  Fig.  13,  is 
for  the  atmospheric  pressure,  ethylene  flame^,  where  the  peak  temperature  is  about  1650K 
and  the  C/O  ratio  is  0.8.  The  analogous  comparison  for  C/O  =  0.96  is  shown  in  Fig.  14. 
While  the  overall  agreement  of  the  sectional  model  is  slightly  better,  the  agreement  of  the 
Lindsfedt  model  is  satisfactory  considering  its  simplicity.  For  the  higher  temperature 
propane  and  acetylene  flames^,  however,  where  peak  temperatures  are  about  1850K  and 
2050K,  respectively,  the  simplified  model  fails,  as  shown  in  Figs.  15  and  16.  All  the 
sectional  model  calculations  were  made  with  the  "modified  Frenklach-Wang"  surface 
growth  mechanism^;  it  incorporates  high  temperature  active  site  decay  and  acetylene 
desorption  mechanisms  that  make  possible  the  better  agreement  shown.  There  is  no  reason 
to  expect  better  performance  from  other  global  models  in  the  literature,  and  it  is  probably 
the  case  that  their  usefulness  is  limited  to  restricted  parameter  ranges  where  they  have  been 
fitted  to  opposed  jet  data,  for  example.  In  these  premixed  flames,  the  global  model  has 
considerable  sensitivity  to  inception  rate,  as  has  been  found  for  the  sectional  model H. 

II.  G.  Coupling  of  Soot  Growth  Model  and  Opposed  Jet  Flow  Code 

The  soot  growth  equations  have  been  programmed  in  a  seamless  fashion  into  an  opposed 
jet  flow  solver,  and  the  combined  program  is  being  debugged.  The  dynamical  equations 
for  the  particulate  concentrations  include  nucleation  coupled  to  the  benzene  production  rate, 
particle  coalescence,  surface  growth,  and  oxidation.  This  work  marks  the  first  time  that  the 
particle  sectional  growth  equations  widely  used  in  aerosol  science  have  been  coupled  to  the 
conservation  equations  for  a  diffusion  flame.  Particle  diffusion  and  thermophoresis 
velocities  in  the  free  molecule  regime  have  been  included.  The  first  set  of  calculations  has 
been  based  on  a  perturbation  approximation  in  which  particle  thermochemistry  and 
depletion  of  gas  phase  species  like  acetylene  by  surface  growth  are  assumed  to  be  small. 
The  next  stage  will  include  these  effects,  as  well  as  radiative  loss.  Preliminary  indications 
are  that  the  addition  of  the  particle  sectional  equations  does  not  complicate  the  convergence 
behavior  of  the  Newton-Raphson  solution  scheme.  When  debugging  is  completed, 
comparisons  with  opposed  jet  data  will  be  made. 
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Fig.  13  Soot  Predictions  for  Premixed 
Ethylene  Flame  (C/0  =  0.8) 


Fig.  14  Soot  Predictions  for  Premixed 
Ethylene  Flame  (C/0  =  0.96) 
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in.  LIST  OF  PUBLICATIONS 


A  paper  entitled  "Influence  of  Radiative  Loss  on  Nitric  Oxide  Formation  in  Counterflow 
Diffusion  Flames  at  High  Pressure",  by  A.  Vranos  and  Robert  J.  Hall,  has  been  published 
in  Combustion  and  Flame.  A  copy  of  this  document,  which  was  partially  supported  under 
this  contract,  is  provided  in  Appendix  B. 

A  paper  entitled  "The  Radiative  Source  Term  for  Plane-Parallel  Layers  of  Reacting 
Combustion  Gases",  by  Robert  J.  Hall,  has  been  published  in  the  Journal  of  Quantitative 
Spectroscopy  and  Radiative  Transfer.  A  copy  of  this  document  is  provided  in  Appendix  D. 

A  manuscript  entitled  "Radiative  Dissipation  in  Gas-Soot  Mixtures",  by  Robert  J.  Hall,  has 
been  submitted  to  the  Journal  of  Quantitative  Spectroscopy  and  Radiative  Transfer.  A  copy 
of  this  document  is  provided  in  Appendix  E. 

A  manuscript  entitled  "The  Formation  of  Benzene  From  Methane  in  Shock  Tubes  and 
Diffusion  Flames"  by  M.  B.  Colket  and  M.  D.  Smooke  is  in  preparation. 

A  review  article  describing  state-of-the-art  modeling  of  soot  formation  processes  is  being 
prepared.  It  is  co-authored  by  I.  Kennedy  and  M.  Colket. 


IV.  MEETING  INTERACTIONS  AND  PRESENTATIONS 

A  paper  entitled  "Predictions  of  Soot  Particle  Growth  Based  on  Aerosol  Dynamics 
Simulations",  by  Robert  J.  Hall  and  Meredith  B.  Colket,  was  presented  at  the  European 
Aerosol  Conference,  September  7-11, 1992,  in  Oxford,  U.  K.  The  extended  abstracts 
were  published  in  a  special  issue  of  the  Journal  of  Aerosol  Science,  and  this  is  presented  in 
Appendix  F. 

Robert  J.  Hall  presented  a  paper  entitled  "Radiative  Transfer  in  Sooting  Counterflow 
Flames"  at  the  Central  and  Eastern  States  Combustion  Institute  Joint  Technical  Meeting, 
March  15-17, 1993,  in  New  Orleans,  LA.  A  copy  of  the  extended  abstract  is  provided  in 
Appendix  C. 

A  paper  entitled  "The  Formation  of  Benzene  From  Methane"  by  M.  B.  Colket  and  M.  D. 
Smooke  was  presented  to  the  Joint  Technical  Meeting  of  the  Centra  States  and  Eastern 
States  Section  of  the  Combustion  Institute,  March  15-17, 1993  in  New  Orleans.  A  copy  of 
the  extended  abstract  is  provided  in  Appendix  A. 

Robert  J.  Hall  made  a  presentation  to  NASA  Lewis  personnel  on  March  25, 1993  that 
included  a  description  of  the  coupled  soot  growth/opposed  jet  code  being  developed  under 
this  contract.  In  a  program  that  will  be  funded  by  NASA,  the  code  will  be  used  to  generate 
sooting  flamelet  libraries  for  a  laminar  flamelet/mixture  fraction  pdf  approach  to  turbulent 
radiation  in  the  rich  zone  of  the  RBQQ  combustor,  under  development  in  the  high  speed 
research  program. 

On  June  2-4, 1993,  M.  Colket  attended  the  15th  Combustion  Research  Conference  held  at 
Lake  Harmony,  Pennsylvania.  M.  Colket  was  invited  to  this  DOE,  Basic  Energy  Sciences 
Contractor’s  Meeting  to  be  an  observer  and  participant. 

V.  RECORD  OF  INVENTIONS 

There  were  no  inventions  during  this  period. 
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The  Formation  of  Benzene  From  Methane 


by  M.  B.  Colket 

United  Technologies  Research  Center 
East  Hartford,  CT  06108 

and 


M.  D.  Smooke 

Department  of  Mechanical  Engineering 
Yale  University 
New  Haven,  CT  06520 


Presentation  to  the  Joint  Technical  Meeting  of  the 
Eastern  States  and  Central  States  Sections  of  the  Combustion  Institute 
March  15-17, 1993,  New  Orleans,  LA 


As  pan  of  an  objective  to  model  soot  production  in  laminar  diffusion  flames,  we  have 
been  analyzing  the  mechanisms  for  the  formation  of  benzene  during  the  oxidative 
pyrolysis  of  methane.  This  presentation  describes  some  experimental  work  on  the  fuel- 
rich  oxidation  of  methane  in  a  single-pulse  shock  tube,  detailed  chemical  kinetic 
modeling  studies  comparing  the  various  routes  to  the  formation  of  benzene,  development 
"of  a  simplified  model  of  benzene  formation  and  destruction  and  finally  model  prediction 
of  benzene  in  opposed-jet  laminar  diffusion  flames. 

A  single-pulse  shock  tube  was  used  to  oxidatively  pyrolyze  methane  over  the  temperature 
range  of  1200  to  2200K,  total  pressures  of  10-12  atmospheres,  equivalence  ratios  of  8  to 
16  and  initial  methane  levels  of  0.6  to  4%  in  argon.  Dwell  times  were  approximately 
500-600  microseconds.  Gas  samples  were  collected  and  analyzed  using  capillary  gas 
chromatography  for  the  hydrocarbons,  hydrocarbon  oxygenates,  carbon  oxides,  oxygen, 
and  hydrogen. 

A  chemical  kinetic  mechanism  for  fuel-rich  methane  oxidation  was  developed  based  on 
the  acetylene  pyrolysis  mechanism  with  added  reactions  for  methane  pyrolysis  and 
oxidation  and  including  reactions  involving  Cs-hydrocarbons  as  well  as  six  separate  steps 
for  the  formation  of  benzene  (or  phenyl).  These  reactions  are: 


A 

F 

Reactions 

(cc/mole/see) 

(Kcal/tuole) 

n-C4H5+C2H2-»C6H6+H 

4.5x10*2 

10.0 

C2H3+C4H4— >C4H6+H 

4.0x10* 1 

0.0 

C2H2+n-C4H3-*l-phenyl+H 

5.2x10** 

0.0 

C3H3+C3H3-»phenyl+H 

5.0x10*2 

0.0 

CH3+c-C5H5->C6H6+2H 

5.0x10*2 

0.0 

C2H2+C-C5H5— »C7H7 

3.0x10** 

15.0 

The  Iasi  step,  which  forms  the  benzyl  radical,  could  lead  to  the  formation  of  toluene  and 
hence  benzene  following  substitution  of  the  methyl  radical.  Evidence  for  the  existence  of 
the  last  two  steps  comes  from  recent  pyrolysis  studies  on  cyclopentadiene  pyrolysis.  The 
reader  should  be  aware  that  despite  the  :  -  lusion  of  the  propargyl  radical  recombination 
step  in  this  study,  recent  results  by  Kiefer  (1992)  place  this  particular  step  in  question  and 
instead  support  the  overall  process  of  propargyl  addition  to  allene  to  explain  the  observed 
facile  formation  of  benzene. 

A  comparison  of  some  of  the  experimental  data  with  predictions  from  the  detailed 
chemical  model  (including  effects  of  the  quenching  wave)  are  shown  in  Fig.  1.  The 
comparisons  aren’t  perfect,  but  at  least  are  reasonable.  A  detailed  reaction  path  analysis 
demonstrated  that  at  1800K  near  the  peak  in  the  benzene  concentration,  the  dominant 
path  for  the  formation  of  aromatics  was  propargyl  recombination.  The  net  rates  of 
reactions  for  each  of  six  aromatic-forming  reactions  is  shown  in  Fig.  2.  The  figure  shows 
the  complexity  in  this  reaction  system.  At  1800K,  although  C3H3  recombination  is 
dominant,  it  is  also  clear  that  decomposition  of  phenyl  (eventually)  into  acetylene  and  n- 
C4H3  is  an  important  removal  process.  Presumably,  at  slightly  lower  temperatures  and 
perhaps  with  different  dominant  reactive  species,  net  formation  of  phenyl  results  from 
this  reaction.  Also  of  interest  from  Fig.  2  is  the  apparent  relative  importance  of  methyl 
radical  addition  to  cyclopentadienyl  radicals.  Obviously  the  proposed  reaction  is  not 
elementary  and  probably  requires  a  reaction  barrier.  The  suggested  rate  constant  is 
consistent  with  an  activation  of  only  5  kcal/mole.  Consequently,  this  reaction  can  be 
considered  speculative. 

Assuming  that  propargyl  recombination  is  in  fact  the  dominant  ring-forming  reaction  for 
the  shock  tube  studies  and  further  assuming  that  this  conclusion  can  be  extended  to 
diffusion  flame  conditions  (See  subsequent  paragraph  for  further  discussion),  then  a 
simplified  kinetic  mechanism  was  constructed.  This  mechanism  was  selected  to  mate 
with  an  existing  Cl  and  C2  mechanism  for  methane  combustion  already  in  use  for 
diffusion  flame  modeling.  This  reduced  set  was  selected  to  minimize  additional  reactions 
and  species.  This  eighteen  reaction  set  includes  eight  new  species  not  present  in  the  C 1-, 
C2  reaction  set,  as  seen  in  Table  I.  Using  this  reaction  set,  benzene  profiles  have  been 
predicted  in  an  opposed  jet  diffusion  flame.  Concentrations  of  a  variety  of  species  are 
shown  in  Fig.  3  for  a  flame  with  mild  stretch.  Peak  benzene  profiles  are  shown  in  Fig.  4 
for  a  range  of  flame  stretch  conditions.  As  the  strain  rate  begins  to  increase,  benzene 
concentrations  rapidly  decrease,  then  there  is  a  more  gradual  decrease  in  the  peak 
benzene  concentration  at  higher  strain  rates.  Qualitatively,  this  result  is  consistent  with 
the  observations  of  Du,  Axelbaum  and  Law,  who  found  the  soot  extinction  limit  to  be 
about  a  factor  of  ten  less  than  that  at  extinction  and  the  fluorescence  limit  to  be  only  a 
factor  of  two  and  one-half  below  the  extinction  limit. 

As  stated  earlier,  the  above  predictions  are  based  upon  the  premise  that  propargyl 
recombination  is  the  dominant  ring-formation  mechanism.  The  modeling  results  indicate 
that  benzene  concentration  peaks  at  temperatures  near  1400K,  significantly  lower  than  the 
temperatures  of  the  shock  tube  experiments  where  peak  benzene  concentrations  were 
observed.  The  temperature  difference  and  the  possible  ’dwell  time'  difference  could  lead 
to  a  difference  in  the  dominant  reaction  pathway.  Consequently,  the  conclusions  reached 
in  this  analysis  should  be  considered  preliminary. 
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Influence  of  Radiative  Loss  on  Nitric  Oxide  Formation 
in  Counterflow  Diffusion  Flames  at  High  Pressure 

A.  VRANOS 

AB  Research  Associates,  South  Windsor,  CT 

and 

R.  J.  HALL 

United  Technologies  Research  Center,  East  Hartford.  CT 

A  theoretical  analysis  is  given  of  the  effect  of  nonluminous  thermal  radiation  on  the  properties  of 
counterflow  diffusion  flames  at  high  pressure.  The  self-consistent  analysis  includes  an  expression  for  gas  band 
radiant  dissipation  in  the  energy  equation  of  a  counterflow  flame  solver.  NO,  formation  rates  and  other 
properties  are  studied  as  a  function  of  strain  rate  for  adiabatic,  optically  thin,  and  optical  thickness -corrected 
flames.  For  adiabatic  flames,  NO,  concentration  and  flame  temperature  increase  continuously  with  decreasing 
strain  rate.  For  radiating  flames,  a  temperature  level  off  is  exhibited  due  to  the  competing  effects  of  heat  loss 
and  extent  of  reaction  as  strain  rate  (inverse  residence  time)  decreases.  For  the  very  lowest  strain  rates,  up  to 
5%  of  the  flame  enthalpy  is  converted  to  radiation,  and  optical  thickness  corrections  start  to  become 
important.  Certain  factors  that  need  to  be  taken  into  account  in  relating  the  isolated  flamelet  results  to 
practical  gas  turbine  combustors  are  also  discussed. 


tive  loss  has  a  pronounced  effect  on  NO  levels 
for  practical  gas  turbine  operating  conditions  is 
shown  to  depend  strongly  on  the  distribution 
of  low  strain  rate  flamclets. 

ANALYTICAL  METHODS 
Stretched  Flame  Model 

The  flow  field  is  computed  using  a  comprehen¬ 
sive  code  developed  by  Smookc  et  al.  [5-7). 
Detailed  kinetics  and  transport  effects  arc  in¬ 
cluded  in  a  quasi-one  dimensional  analysis  ol 
diffusive  combustion,  namely  air  and  luel  fed 
oppositely  to  the  reaction  /one.  The  structure 
of  the  flame  is  obtained  as  a  solution  ol  a  set 
of  coupled  nonlinear,  two  point,  boundary  value 
problems  along  ihe  stagnation  streamline.  The 
strain  rate  enters  the  calculation  as  a  boundary 
condition  that  specifics  the  velocity  compo¬ 
nents  along  the  edge  of  (he  boundary  layer. 
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INTRODUCTION 

The  objective  of  this  study  was  to  estimate  the 
influence  of  nonluminous  radiative  heat  loss 
on  nitric  oxide  formation  in  high-pressure, 
counterflow  diffusion  flames.  Previous  studies 
of  atmospheric  counterflow  diffusion  flames 
have  shown  that  radiative  cooling  can  reduce 
NO,  concentrations  appreciably  [1-4],  In  this 
report,  an  analysis  of  radiative  heat  loss  in  a 
high  pressure  (10.5  bar)  counterflow  flame  is 
described.  The  resulting  formulation  is  used  as 
a  basis  for  discussing  the  influence  of  radiative 
loss  on  nitric  oxide  formation  in  an  aircraft  gas 
turbine  combustor.  An  order  of  magnitude  es¬ 
timate  of  combustor  mean  strain  rate,  com¬ 
puted  from  one-dimensional  scalar  decay  rate 
measurements  in  a  cross  flow  jet  mixing  con¬ 
figuration,  is  used  to  establish  an  appropriate 
range  of  strain  rates  for  the  counterflow  calcu¬ 
lation.  The  question  whether  gas  band  radia- 
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The  kinetics  data  set,  limited  to  Cl  chemistiy 
for  methane  fuel,  consisted  of  126  reactions 
and  30  species  culled  from  a  more  comprehen¬ 
sive  set  of  213  reactions  and  57  species  devel¬ 
oped  by  Glarborg  et  al.  [8],  N02  chemistry  was 
not  considered.  The  reaction  set  and  kinetic 
rate  constants  are  present  in  Appendix  A.  The 
methane  fuel  temperature  was  selected  to  pro¬ 
vide  the  same  adiabatic  flame  temperature  as 
jet  fuel  at  the  same  operating  conditions.  Com¬ 
bustor  pressure  was  10.5  atm,  typical  of  the 
cruise  condition. 


Theory  of  Radiative  Power  Loss  from 
Counterflow  Flames 


The  effect  of  gas  band  radiation  on  flame 
temperature  and  species  concentrations  was 
calculated  in  a  self-consistent  manner  by  incor¬ 
porating  a  nonadiabatic  radiative  power  loss 
term  for  a  plane-parallel  layer  into  the  energy 
equation  of  the  counterflow  flame  code  [6,  7]. 
In  the  notation  of  Refs.  6  and  7  the  energy 
equation  becomes 


d  l  dT\  dT 
dy  \  dy  )  Cp^  dy 


*  dT 

LpYkVkycpk  — 


*  dqr 

I  *kWkhk  -~=o, 

k-i  dy 


(i) 


where  the  term  dqr/dy  is  the  divergence  of  the 
net  radiative  flux,  qr.  The  nomenclature  for  the 
rest  of  Eq.  1  is  as  given  in  the  references;  the 
other  terms  on  the  l.h.s.  reading  from  left  to 
right  represent  conduction,  convection,  diffu¬ 
sive  transport,  and  energy  release,  respectively. 
The  flux  divergence  or  radiative  source  term 
was  calculated  in  the  manner  given  in  Ref.  9, 
which  used  certain  of  the  results  given  in  Ref. 
10.  It  is  based  on  the  use  of  wideband  gas 
radiation  models  for  H;0.  C02,  and  CO  with 
provision  for  finite  optical  thickness  effects 
in  the  high  pressure  broadening  limit.  Self¬ 
absorption  of  emitted  radiation  was  calculated 
by  taking  hemispherical  averages  of  the  solu¬ 
tion  for  the  radiative  flux,  and  integrating  the 
result  over  the  molecular  absorption  band- 
shapes.  Although  the  optically  thin  approxima¬ 


tion  will  usually  be  adequate,  finite  thickness 
effects  can  become  important  at  low  strain  rate 
and  high  pressure,  as  will  be  seen.  Previous 
theoretical  analyses  [1-4]  have  been  limited  to 
the  optically  thin  case.  The  radiative  source 
term  is  accordingly 


ap 

Jr-  ( L(  1  +  oj/w) 
1  +  <o/w 


X  y*(  1  4-  a)/ o>,  r )  +  £j(r) 


O.y-y/ 


xfdy'Ibij  (2) 

/  dy'lbij  =  f  dy'Iblj  (y  >  yf), 

J  J  —  as 

=  /  dy'hi,  (y<>y). 

y 

where  the  first  and  second  terms  on  the  r.h.s. 
represent  emission  and  absorption,  respec¬ 
tively;  the  summation  ij  is  over  all  active  bands 
of  H20,  C02,  and  CO;  represents  the 
integrated  band  intensity;  p,  the  mass  density 
of  the  active  molecule:  fhl]  the  Planck  function 
evaluated  at  the  band  center  frequency;  r  is 
the  band  center  optical  depth  between  the 
point  of  maximum  flame  temperature,  yf,  and 
y;  and  w  is  the  bandwidth  parameter.  Over- 
bars  in  the  nonlocal  absorption  term  denote 
cither  path-averaged  properties  or  properties 
based  on  path-averaged  temperature  taken  be¬ 
tween  the  point  y  and  y, .  I',  and  y*  are 

the  exponential  integral,  gamma  function,  and 
incomplete  gamma  function,  respectively  111]. 
Band  overlap  contributions  are  ignored  in  this 
analysis.  A  quantity  of  interest  that  can  be 
calculated  from  the  converged  solution  is  the 
fraction  of  the  enthalpy  that  is  converted  to 
radiation.  Given  that  the  fourth  term  in  Eq.  1 
is  the  local  rate  of  eneigy  release,  the  Irae- 
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tional  radiative  power  loss  is  [12] 


dqr 

r-'  +  * 


-  If  wk Wkhk  dy | 


(3) 


The  absorption  term  in  Eq.  2  can  represent 
a  computational  burden,  and  the  overall  econ¬ 
omy  of  the  program  was  maintained  by  omit¬ 
ting  the  radiative  terms  from  the  recalculation 
of  the  Jacobian  matrix  employed  in  the  New¬ 
ton’s  method  solution  of  the  nonlinear  dif¬ 
ference  equation  set  [12].  Eq.  2  itself  is  an 
“analytic  approximation”  to  an  exact  but  time- 
consuming  solution  for  the  thickness  term  that 
has  been  suggested  by  the  results  of  Ref.  10.  In 
the  general  formulation  of  the  radiative  source 
term  given  in  Ref.  9,  it  is  possible  to  take  into 
account  nonvacuum  boundary  conditions 
(emission  from  black,  hot  boundaries),  but  this 
was  omitted  from  the  present  calculations. 


Combustor  Mean  Strain  Rate 

The  preceding  analysis  is  used  to  assess  the 
influence  of  gas  band  radiation  on  nitric  oxide 
formation  in  an  aircraft  gas  turbine  combustor. 
The  combustor  flame  is  viewed  as  an  ensemble 
of  flamelets  converted  by  the  mean  flow,  and 
the  combustion  process  is  assumed  to  be  dif¬ 
fusion  driven  (no  premixing)  with  air  and  fuel 
entering  through  separate  boundaries.  Rele¬ 
vant  flamelet  properties  such  as  maximum 
flame  temperature  and  NO  concentration  are 
described  conveniently  as  a  function  of  the 
strain  rate  boundary  condition.  For  otherwise 
fixed  boundary  conditions,  the  flamelet  strain 
rate  and  mean  flamelet  scalar  dissipation  rate 
are  uniquely  related.  Thus,  the  flamelet  model 
can  be  applied  to  a  combustor  flame  if  the 
mean  combustor  scalar  dissipation  is  known. 
An  order  of  magnitude  estimate  of  the  dissipa¬ 
tion  rate  is  obtained  from  recent  measure¬ 
ments  of  unmixedness  in  a  cylindrical  model 
combustor  with  multijet  cross  flow  injection 
[13].  Planar  distributions  of  unmixedness  are 
used  to  compute  area-average  unmixedness  as 
a  function  of  the  streamwise  coordinate.  The 


one-dimensional,  streamwise  rate  of  decay  of 
jet  fluid  concentration  is  inferred  from  the 
mean  velocity. 

The  following  equations  describe  the  volu¬ 
metric  mean  scalar  dissipation  rate  for  normal 
injection  from  a  row  of  holes  in  a  single  plane 
normal  to  the  duct  axis. 


X(X/D)  =  0.52 


—J 


(1  ~Z) 


-0.654 


U/D 


(X/D)  o  =  o.4or0  203. 


where  \  isjhe  volumetric  mean  scalar  dissipa¬ 
tion  rate,  Z  is  the  mass  fraction  of  jet  fluid 
based  on  the  jet  and  approach  stream  mass 
flows,  U  is  the  approach  stream  mean  velocity, 
and  J  is  the  single  jet  momentum  flux  ratio. 

It  is  seen  that  x  depends  on  Z,  and  so  is,  in 
fact,  a  function  of  the  number  of  injection 
stages  and  the  rate  of  air  addition  per  stage.  In 
this  example,  it  is  assumed  that  one  half  the 
dilution  air  is  added  at  the  primary  zone  exit 
and  that  this  fluid  is  thoroughly  mixed  in  one 
combustor  height  (or  diameter).  For  typical 
values  of  J  =  23  and  Z  =  0.53,  x  =  121  s'1 
corresponding  to  a  mean  strain  rate  of  330  s  1 . 

This  estimate  of  mean  strain  rate  can  be 
checked  for  correct  order  of  magnitude  using 
the  fundamental  relationship  for  viscous  dissi¬ 
pation  [14], 


e„  =  A(u')3/ A, 


(5) 


where  A  is  a  constant  near  unity  (herein  as¬ 
sumed  to  be  unity),  //'  is  the  rms  turbulence 
intensity  for  isotropic  turbulence,  and  A  is  the 
integral  scale  of  turbulence.  Assuming  typical 
values  of  u’  =  129  cm/s  and  A  -  1.0  cm  for 
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the  combustor  primary  zone,  the  resulting 
flameiet  dissipation  is  2.1  x  10 6  cm2/s\  A 
mean  dissipation  within  the  flameiet  is  given  by 
the  expression 


where  v  is  the  average  flameiet  kinematic  vis¬ 
cosity.  Substituting  v0  =  -lay,  where  a  is  the 
strain  rate  and  v0  is  the  normal  velocity  rela¬ 
tive  to  the  stagnation  point,  and  letting  ev  =  ef 
yields  a  =  474  s'1. 

RESULTS 

The  influence  of  radiative  loss  on  temperature 
is  seen  in  Fig.  1,  where  peak  flameiet  tem¬ 
perature  is  plotted  as  a  function  of  strain  rate 
for  adiabatic,  optically  thin  and  thickness- 
corrected  calculations.  The  effect  of  radiation 


on  temperature  becomes  more  pronounced 
as  strain  rate  decreases,  where  it  can  also 
be  seen  that  thickness  corrections  have  to  be 
taken  into  account  for  the  lowest  strain  rates 
(20-30  s'1).  All  curves  converge  at  sufficiently 
high  strain  rate  as  a  result  of  diminishing  ra¬ 
diative  loss.  These  results  are  mirrored  in  the 
dependence  of  peak  nitric  oxide  concentration 
on  strain  rate,  shown  in  Fig.  2.  The  sensitivity 
to  temperature  is  in  accordance  with  the  fact 
that  most  of  the  nitric  oxide  is  formed  via  the 
thermal  (Zeldovich)  mechanism.  Differences 
between  the  optically  thin  and  thickness-cor¬ 
rected  model  also  become  more  apparent  at 
low  strain  rate.  The  corresponding  radiative 
loss,  expressed  as  a  fraction  of  the  integrated 
energy  release  is  presented  in  Fig.  3.  Within 
the  framework  of  the  present  analysis,  it  is  not 
possible  to  make  a  definite  statement  that  non- 
luminous  radiation  should  have  a  significant 
influence  on  nitric  oxide  formation  in  a  practi- 


Fig.  I  Peak  temperature  in  high- pressure  cmimerfli m  flame,  icrvi.  .'i.iiii  rale  with  and  without  radiation  The  range 
of  strain  rales  is  20-KIHl  s  1 
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Log(lO)  Strain  rate  -  sec  1 

Fig.  2.  Effect  of  radiation  on  peak  NO  concentrations  in  high-pressure  counterflow  flames. 


cal  application.  However,  if  the  scalar  dissipa¬ 
tion  rate  is  log-normally  distributed,  the  calcu¬ 
lations  suggest  that  significant  contributions  to 
the  weighted  sum  of  nitric  oxide  could  be 
expected  from  flamelets  whose  scalar  dissipa¬ 
tion  rates  are  below  the  mean  value.  For  exam¬ 
ple,  at  5%  loss,  the  reduction  in  peak  NO 
levels  relative  to  the  adiabatic  solution  is  about 
30%.  Thus,  it  appears  that  detailed  calcula¬ 
tions  of  radiation  in  combustors  are  worth 
pursuing.  In  this  regard,  practical  applications 
require  consideration  of  the  influence  of 
neighboring  flamelets  on  net  emission  from  the 
turbulent  flame.  Radiative  interactions  be¬ 
tween  flamelets  and  overall  optical  thickness 
effects  could  cause  gas  cooling  rates  in  practi¬ 
cal  contexts  to  be  less  than  that  calculated 
here.  Radiative  loss  from  soot  and  its  effect  on 
nitric  oxide  could  be  much  more  significant. 
The  theoretical  expression  for  the  radiative 
loss  (Eq.  2)  can  be  extended  readily  to  include 
particulate  emission.  There  is  a  need,  too,  for  a 
more  comprehensive  and  systematic  study  of 


the  effects  of  pressure,  peak  temperature,  and 
strain  rate  on  radiative  loss  in  opposed  jet 
flames.  This  is  beyond  the  scope  of  the  present 
article,  however,  and  will  be  addressed  in  a 
future  publication. 

CONCLUSIONS 

A  gas  band  radiation  term  has  been  included 
in  the  energy  equation  for  counterflow  diffu¬ 
sion  flames,  providing  a  self-consistent  analysis 
of  the  effects  of  radiation  on  flame  tempera¬ 
ture  and  NO  levels.  For  conditions  representa¬ 
tive  of  a  gas  turbine  combustor,  radiative  cor¬ 
rections  can  become  important  for  low  strain 
rates,  with  up  to  5%  conversion  of  enthalpy  to 
radiation,  and  reductions  of  approximately  30rv 
in  peak  NO  levels  relative  to  the  adiabatic 
solution.  At  high  pressures  and  low  flamelct 
strain  rates,  optical  thickness  corrections  arc 
required  for  high  accuracy.  The  question 
whether  gas  band  radiative  loss  has  a  pro¬ 
nounced  effect  on  NO  levels  for  practical  gas 
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Fig.  3.  Nonadiabatic  radiative  power  loss  in  counterfiow  flames  versus  strain  rate. 


turbine  operating  conditions  is  shown  to  de¬ 
pend  strongly  on  the  distribution  of  low- 
strain-rate  flamelets. 
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APPENDIX  A 


Reaction  Mechanism  Rate  Coefficients  in  the  Form  kf  -  AT®  exp(  -£0  /  RT ).  Units  are  moles,  cubic  centimeters, 

seconds,  Kelvin,  and  calories  /  mole. 


Reaction 

A 

P 

E 

1. 

H  +  02  -  O  +  OH 

5.10E16 

-0.820 

16510. 

2. 

H2  +  O  -  H  +  OH 

1.80E10 

1.0 

8830. 

3. 

H2  +  OH  =  H20  +  H 

1.20E09 

1.3 

3630. 

4. 

OH  +  OH  =  H20  +  O 

6.0E08 

1.3 

0. 

5. 

H  +  OH  +  M  =  HjO  +  M  (M  =  AR) 

H2O/20./ 

7.50E23 

-2.6 

0. 

6. 

o2  +  m  =  o  +  o  +  m 

1.90E11 

0.5 

95560. 

7. 

H  +  H  +  M  =  H2  +  M  (M  =  AR) 

h2o/o.o/h2/o.o/co2/o.o/ 

1.0E18 

-1.0 

0. 

8. 

H  +  H  +  H2  =  H2  +  H2 

9.20E16 

-0.6 

0. 

9. 

H  +  H  +  H20  =  H2  +  HjO 

6.00E19 

-1.250 

0. 

10. 

H  +  H  +  C02  =  H2  +  C02 

5.49E20 

-2.0 

0. 

11. 

H2  +  02  =  OH  +  OH 

1.70E13 

0.0 

47780. 

12. 

H  +  02  +  M  «  H02  +  M  (M  -  AR) 
H2O/21./COj/5./H2/3.3/CO/2./O2/0./N2/0./ 

2.10E18 

-1.0 

0. 

13. 

H  +  02  +  02  =  H02  +  02 

6.70E19 

-1.420 

0. 

14. 

H  +  02  +  N2  -  H02  +  N2 

6.70E19 

-  1.420 

0. 

15. 

ho2  +  h  =  h2  +  o2 

2.50E13 

0. 

700. 

16. 

H02  +  H  -  OH  +  OH 

2.50E14 

0. 

1900. 

17. 

H02  +  O  =  OH  +  02 

4.80E13 

0. 

1000. 

18. 

H02  +  OH  -  H20  +  02 

5.00E13 

0. 

1000. 

19. 

ho2  +  ho2  -  h2o2  +  o2 

2.00E12 

0. 

0. 

20. 

h2o2  +  m  -  oh  +  oh  +  m 

1.20E17 

0. 

45500. 

21. 

h2o2  +  H  =  ho2  +  h2 

1.70E12 

0.0 

3750. 

22. 

H2Oj  +  OH  -  HjO  +  HOj 

1.00E13 

0.0 

1800. 

23. 

CO  +  O  +  M  -  C02  +  M 

3.20E13 

0.0 

-4200. 

24. 

CO  +  Oj  *  co2  +  o 

2.50E12 

0. 

47700. 

25. 

CO  +  OH  *  C02  +  H 

1.50E07 

1.3 

-760. 

26. 

CO  +  H02  -  C02  +  OH 

5.80E13 

0. 

22930. 

27. 

CH4  +  M  -  CH,  +  H  +  M  (M  -  AR) 

H20/5./ 

1.00E17 

0. 

88000. 

28. 

CH4  +  H  =  CH,  +  H2 

2.20E04 

3.0 

8750. 

29. 

CH4  +  O  =  CH,  +  OH 

1.20E07 

2.080 

7630. 

30. 

CH4  +  OH  =  CH  3  +  HjO 

3.50E03 

3.080 

2000. 

31. 

CH4  +  CH,  =  CH,  +  CH, 

1.30E13 

0. 

9500. 

32. 

CH,  +  M  -  CH2  +  H  +  M 

1.90E16 

0. 

91600. 

33. 

CH,  +  H  =  CH2  +  H2 

9.00E13 

0. 

15100. 

34. 

CH,  +  O  -  CH20  +  H 

6.80E13 

0 

0. 

35. 

CH,  +  O  =  CH,  +  OH 

5.00E13 

0. 

12000. 

36. 

CH,  +  OH  =  CH,  +  H20 

1  JOE  13 

0. 

5000. 

37. 

CH,  +  OH  =  CH  ,0  +  H2 

1.00EI2 

0 

0 

38. 

CH,  +  O,  =  CH,0  +  OH 

5.20E13 

0. 

34570 

39. 

CH,  +  O,  =  CH  ,0  +  O 

7.00E12 

0 

25650. 

40. 

CH  ,0  +  M  =  CH  ,0  +  H  +  M 

I.00EI4 

0 

25000 

41. 

CH  ,0  +  H  =  CH20  +  H2 

2.00EI3 

0.0 

0. 

42. 

CH,0  +  O  =  CH,0  +  OH 

1.00E13 

0 

(V 

43. 

CH  ,0  +  OH  =  CH,0  +  H,0 

1.00E13 

0. 

0. 

44. 

CH  ,0  +  O,  =  CH,0  +  HO, 

6.30EI0 

0 

2600. 

45. 

CH,0  +  M  =  HCO  +  H  +  M 

3.31E16 

0.0 

81000. 

46. 

CH,0  +  H  =  HCO  +  lf2 

2.20EII8 

1.77(1 

10500 

47. 

CH  ,0  +  O  =  HCO  +  OH 

1  801  13 

0.0 

30X0 

48. 

CH  ,0  +  OH  =  HCO  +  H,0 

3  401(C) 

1  180 

447. 

49. 

HCO  +  M-CO  +  H  +  M 

1.60E14 

0 

14700. 
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50. 

HCO  +  H  =  CO  +  H; 

4.00E13 

0.0 

0.0 

51. 

HCO  +  O  =  CO  +  OH 

3.00E13 

0.0 

0.0 

52. 

HCO  +  O  =  C02  +  H 

3.00E13 

0.0 

0.0 

53. 

HCO  +  OH  =  CO  +  HjO 

5.00E12 

0.0 

0.0 

54. 

HCO  +  O;  =  CO  +  H02 

3.30E13 

-0.4 

0.0 

55. 

ch2  +  H  =  CH  +  h2 

7.30E17 

-1.560 

0. 

56. 

ch2  +  o  =  co  +  h  +  h 

3.00E13 

0.0 

0. 

57. 

CH,  +  O  =  CO  +  H2 

5.00E13 

0.0 

0. 

58. 

CH2  +  O  =  CH  +  OH 

5.00E13 

0.0 

12000. 

5° 

CH2  +  OH  =  CHjO  +  H 

3.00E13 

0.0 

0. 

60. 

CH2  +  OH  =  CH  +  HjO 

4.50E13 

0.0 

3000. 

61. 

CH2  +  02  =  C02  +  H  +  H 

1.60E12 

0.0 

1000. 

62. 

CH2  4-  02  =  C02  +  H2 

6.90E11 

0.0 

500. 

63. 

ch2  +  o2  =  CO  +  h2o 

1.90E10 

0.0 

-1000. 

64. 

ch2  +  o2  =  CO  +  oh  +  h 

8.60E10 

0.0 

-500. 

65. 

CH,  +  02  =  HCO  +  OH 

4.30E10 

0.0 

-500. 

66. 

CH]  +  02  =  CH,0  4-  O 

2.00E13 

0.0 

9000. 

67. 

ch2  +  co2  -  co  +  ch2o 

1.10E11 

0.0 

1000. 

68. 

CH  +  O  =  CO  +  H 

5.70E13 

0.0 

0. 

69. 

CH  +  OH  *=  HCO  +  H 

3.00E13 

0.0 

0. 

70. 

CH  +  O,  =  HCO  +  O 

3.30E13 

0.0 

0. 

71. 

CH  +  CO,  =  HCO  +  CO 

3.40E12 

0.0 

690. 

72. 

NHj  +  M*-  NH2  +  H  +  M 

1.40E16 

0.0 

90600. 

73. 

NH3  +  H  -  NH2  +  H2 

7.00E06 

2.390 

10171. 

74. 

NH,  +  O  =  NH,  +  OH 

2.10E13 

0.0 

9000. 

75. 

NH,  +  OH  =  NH2  +  H20 

2.04E06 

2.04 

566. 

76. 

NH2  +  H  «  NH  +  H2 

6.90E13 

0.0 

3650. 

77. 

NH2  +  O  -  NH  +  OH 

6.80E12 

0.0 

0.0 

78. 

NH2  +  O  -  HNO  +  H 

6.60E14 

-0.5 

0. 

79. 

NH2  +  OH  =  NH  +  H20 

4.50E12 

0.0 

2200. 

80. 

NH,  +  N-  Nj  +  H  +  H 

7.20E13 

0.0 

0 

81. 

NH*,  +  NO  -  N,  +■  H  +  OH 

8.80E15 

-1.250 

0. 

82. 

NH2  +  NO  =  N2  +  H20 

3.80E15 

-1.25 

0.0 

83. 

NH  +  H  -  N  +  H, 

3.00E13 

0.0 

0. 

84. 

NH  +  O  =  NO  +  H 

2.00E13 

0.0 

0. 

85. 

NH  +  OH  =  HNO  +  H 

2.00E13 

0.0 

0. 

86. 

NH  +  OH  =  N  +  H,0 

5.00E11 

0.5 

2000. 

87. 

NH  +  O,  =  HNO  +  0 

1.00E13 

0.0 

12000. 

88. 

NH  +  02  =  NO  +  OH 

I.4E1 1 

0.0 

2000. 

89. 

NH  +  NO  =  N,0  +  H 

4.30E14 

-0.5 

0.0 

90. 

NH  *  N  =  N.  +  H 

3.00E13 

0.0 

0. 

91. 

N  +  O,  =  NO  +  O 

6.40E09 

1.0 

6280. 

92. 

N  +  OH  =  NO  +  H 

3.80E13 

0.0 

0.0 

93. 

N  +  NO  =  N.  +  O 

3.30E12 

0.3 

0. 

94. 

N  +  CO.  =  NO  +  CO 

1.90E11 

0.0 

3400. 

95. 

HNO  +  M  =  H  +  NO  +  M  (M  =  AR) 

H.O/6  /II./2  /0./2./N,/2./ 

1.50E16 

0.0 

48680. 

96 

UNO  *  II  11  -  NO 

5.00E+  12 

0.0 

0.0 

97. 

HNO  ♦  OH  NO  *  H.O 

3.6E13 

0.0 

0.0 

98. 

N.O  •  M  N-  t  O  +  M 

I.60E14 

0.0 

51600 

99. 

N,0  -  II  -  N  *  OH 

7.60E13 

0.0 

15200. 

100 

N.O  *  O  -  NO  +  NO 

1.00E14 

0.0 

28200. 

101 

N.O  *  O  *  N.  *  0; 

I  00E14 

0.0 

28200. 

102. 

HCN  *  O  =  CN  4  OH 

2.70E09 

1.580 

26600 

103. 

HCN  *  C)  -  NCO  4  II 

1  40E04 

2.64 

4980. 

104 

HCN  *  O  -  NH  *  CO 

3.50E03 

2.64 

4980 

105 

HCN  *  OH  CN  *  11,0 

1.5E13 

0.0 

10929. 

106 

CN  -  O  CO  ,  N 

I.80E13 

0.0 

0 

107 

CN  -  OH  VO  *  H 

6.(K)E13 

0.0 

0.0 

108 

CN  *11.4  HCN  4  H 

3.00E05 

2.45 

2237. 

109. 

CN  4  0.  -  NCO  4  O 

5.60EI2 

0.0 

0.0 
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no. 

CN  +  NjO  =  NCO  +  N, 

1.00E13 

0.0 

0. 

111. 

NCO  +  M  =  N  +  CO  +‘m 

3.10E16 

-0.5 

48000. 

112 

NCO  +  H  =  NH  +  CO 

5.00E13 

0.0 

0. 

113. 

NCO  +  O  =  NO  +  CO 

5.60E13 

0.0 

0. 

114. 

NCO  +  OH  =  NO  +  CO  +  H 

1.00E13 

0.0 

0. 

115. 

NCO  +  N  =  N2  +  CO 

2.00E13 

0.0 

0. 

116. 

NCO  +  NO  =  N20  +  CO 

1.00E13 

0.0 

-390. 

117. 

CH  +  NO  =  HCN  +  O 

1.10E14 

0.0 

0. 

118. 

CH  +  N2  =  HCN  +  N 

1.90E11 

0.0 

13600. 

119. 

CH2  +  N2  =  HCN  +  NH 

1.00E13 

0.0 

74000. 

120. 

CH  +  NH2  =  HCN  +  H  +  H 

3.00E13 

0.0 

0. 

121. 

CH  +  NH  =  HCN  +  H 

5.00E13 

0.0 

0. 

122. 

CH2  +  NH  -  HCN  +  H  +  H 

3.00E13 

0.0 

0. 

123. 

CH  +  N  -  CN  +  H 

1.30E13 

0.0 

0. 

124. 

CH2  +  N  =  HCN  +  H 

5.00E13 

0.0 

0. 

125. 

CH,  +  N  =  HCN  +  H  +  H 

5.00E13 

0.0 

0. 

126. 

CH4  +  N  =  NH  +  CH3 

1.00E13 

0.0 

24000. 

“Third-body  efficiencies:  ks( H20)  —  20/fcs(Ar). 

6 Third-body  efficiencies:  &7(H20)  =  /c,(H2)  =  lc7(C02)  =  0<:7(Ar). 

“Third-body  efficiencies:  A:12(H20)  =  2U12(Ar),  *12(C02)  =  5fc12(Ar),  Jt,2(H2)  =  3.3ik12(Ar),  lcI2(CO)  =  2Jt12(Ar), 
*12(02)  -  0*I2(Ax),  *12(N2)  -  0*,2(Ar). 

‘'Third-body  efficiencies:  Jk27(H20)  «  5/fc27(Ar). 

'Third-body  efficiencies:  i^(H20)  =  6*,5(Ar),  /c95(H2)  =  k9!(02)  =  *,5(N2)  =  2k9S(Ai). 
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Radiative  Transfer  in  Sooting  Counterflow  Flames 

by 


Robert  J.  Hall 

United  Technologies  Research  Center 
East  Hartford,  Conn. 


It  is  well  known  that  substantial  fractions  of  flame  energy  can  be  converted  to  radiation  (Ref. 
1),  and  that  the  gas  cooling  resulting  from  this  non-adiabatic  loss  can  in  turn  affect  flame  chemistry, 
and  lead  to  inaccuracies  in  prediction  of  pollutants  like  nitric  oxide.  These  effects  are  expected  to 
be  most  extreme  in  sooting  flames.  Quantitative  assessment  of  these  effects  requires  inclusion  of 
an  energy  sink  term  in  the  flow  energy  equation.  This  sink  term  is  given  by  the  divergence  of  the 
net  radiative  flux,  which  is  the  difference  of  an  emission  term  valid  in  the  optically  thin  limit,  and 
a  self-absorption  term  that  becomes  important  when  optical  thickness  effects  need  to  be  accounted 
for.  The  expression  for  the  emission  term  is  generally  trivial  even  when  both  soot  and  gas  band 
radiation  are  present,  and  does  not  depend  on  flame  geometry  because  only  local  properties  are 
required.  The  self-absorption  term  is  generally  more  difficult  because  a  solution  for  the  radiative 
flux  is  implicitly  required.  This  term  is  complicated  further  when  both  soot  and  gas  radiators  are 
present  because  of  absorption  profile  overlaps.  In  this  paper,  a  complete  solution  for  the  radiative 
dissipation  is  given  for  the  one-dimensional,  plane-parallel  or  boundary  layer  flow  problem.  Both 
gas  band  and  soot  radiation  are  considered,  and  the  resulting  solution  is  valid  for  all  degrees  6f 
optical  thickness.  The  expression  is  given  in  semi-analytic  form  in  terms  of  single  quadratures 
normal  to  the  flame  structure.  The  theory  will  be  illustrated  by  example  calculations  for  a  model, 
high  pressure,  sooting  counterflow  diffusion  flame. 


The  main  features  of  the  theoretical  analysis  are  as  follows.  The  absorption  terms  are  derived  by 
taking  hemispherical  averages  of  frequency-dependent  solutions  of  the  equation  of  radiative  transfer, 
and  integrating  them  over  the  absorption  bandshapes  of  both  the  molecules  and  the  soot.  The  soot 
absorption  coefficient  is  assumed  to  have  the  Rayleigh  form,  with  scattering  ignored.  For  nonuniform 
paths,  the  soot  terms  can  be  derived  without  approximation;  for  the  g as  band  terms  a  form  of  the 
mean  properties  assumption  is  employed.  The  gas  absorption  profiles  are  assumed  to  have  the 
wideband  model  form,  and  the  transmissivities  are  appropriate  to  the  high  pressure-broadening 
limit.  The  absorptive  part  is  expressible  as  four  terms,  which  are  too  lengthy  to  reproduce  here. 
Only  the  pure  soot  solution  will  thus  be  given  here.  Assuming  that  y  represents  the  coordinate 
normal  to  the  flame  structure,  the  emission  and  absorption  terms  of  the  net  radiative  dissipation  Q 
=  Qemiss  -  Qabs  are  expressible  as: 


coo  4 

Qemiss  =  4*c,  5ci  ■  /  r  du  •  cf  fv(y)T5(y) 

Jo  e  -  1 


(1) 


where  ci  and  ct  are  the  Planck  function  constants,  c*  is  the  constant  in  the  Rayleigh  soot  absorption 
coefficient  law  K,b.(u>)  =  c,wfv ,  fv  is  the  local  soot  volume  fraction,  T  is  the  local  temperature,  and 
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Qabs  =  2>r  cl  c,  f,(y) 


5! 


dy'  My')  T°(y') 


£» 


where 


P..(y'h 

Q»(y')' 


V  1  (  Pi>(y,) 

^  k  Vpn(y')  +  Q„(y')/ 


P. .(y')  =  nc2/T(y') 

Q. .(y')  =  c,fv  |y  -  y'l 


(2) 


An  implicit  assumption  lias  been  that  the  dispersion  of  the  soot  index  of  refraction  is  small. 
Other  terms  give  the  rate  of  gas  absorption  of  gas  radiation,  gas  absorption  of  soot  radiation,  and 
soot  absorption  of  gas  radiation,  and  are  expressible  in  terms  of  exponential  integrals. 

The  generalized  net  radiative  dissipation  for  both  soot  and  gas  molecules  has  been  incorporated 
into  a  counterflow  diffusion  code  solver  ,  which  includes  detailed  kinetics  and  transport  effects  (Ref. 
2).  The  concentrations  of  the  radiating  species  COj,  HjO,  and  CO  are  derived  in  the  counterflow 
solution,  but  the  soot  has  been  treated  parametrically.  Example  calculations  have  been  performed 
for  a  10.5  atm,  methane-air  flame.  The  pressure  approximates  that  of  an  aircraft  gas  turbine 
combustor  at  the  cruise  condition,  and  the  fuel  temperature  has  been  adjusted  to  achieve  a  peak 
temperature  approximating  the  adiabatic  flame  temperature  of  jet  fuel  (Ref.  3).  The  adiabatic  (no 
radiation)  temperature  profile  for  this  flame  is  shown  in  Figure  1,  along  with  the  assumed  location 
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of  the  soot  profile.  Since  a  soot  growth  model  has  not  yet  been  included  in  the  counterflow  solver, 
the  location,  peak  soot  volume  fraction,  and  width  of  the  sooting  7,one  are  treated  as  parameters 
uncoupled  from  the  flame  solution,  in  the  calculations  to  follow,  the  peak  in  the  soot  profile  is 
assumed  to  occur  at  about  1650  K,  with  a  bandwidth  of  200  K,  and  the  peak  volume  fraction  is 
varied.  The  strain  rate  in  these  calculations  is  56  sec  *,  a  value  which  is  felt  to  lie  toward  the  lower 
end  of  the  range  of  values  of  practical  interest  (Itcf.  3);  the  strain  rate  enters  the  calculation  as  a 
boundary  condition  which  specifies  the  velocity  components  along  the  edge  of  the  boundary  layer 

Figure  2  shows  the  fractional  (%)  conversion  of  flame  energy  to  radiation  as  the  peak  soot 
volume  fraction  is  varied.  The  lower  asymptotic  limit  represents  the  contribution  of  gas  band 
radiation  alone.  For  peak  soot  volume  fractions  much  in  excess  of  10-c,  the  soot  contributions 
become  important,  and  for  values  approaching  IO-4,  more  than  10%  of  the  flame  energy  is  converted 
to  radiation.  A  peak  soot  volume  fraction  of  this  magnitude  might  not  be  unreasonable  at  high 
pressure.  Losses  of  this  magnitude  cool  the  gas  sufficiently  that  marked  effects  are  seen  on  the 
predicted  nitric  oxide  profiles,  as  shown  in  Fig.  3.  The  kinetics  data  set  is  limited  to  Cl  chemistry 
for  methane  fuel,  and  consists  of  126  reactions  and  30  species  derived  from  a  more  comprehensive 
mechanism  developed  by  Glarborg,  et.  al.  (Ref.  4).  Given  the  strong  temperature  dependences 
exhibited  in  Equations  1  and  2,  the  results  will  clearly  have  considerable  sensitivity  to  the  location 
of  the  soot  line. 


Log(lO)  Peak  Soot  Volume  Fraction 


Figure  2.  Variation  of  percentage  conversion  of  flame  energy  to  radiation  with  peak  soot 
volume  fraction. 
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Figure  3.  Effect  of  soot  loading  on  nitric  oxide  profiles.  The  values  of  volume  fraction  arc 
1,  2,  Si  5x10'°,  1,  2,  Si  5x10-*,  10~4. 

The  full  gas  band  and  soot  radiation  expressions  will  be  given  in  the  paper.  Other  topics  which 
will  be  discussed  include  the  incorporation  of  realistic  radiative  wall  boundary  conditions,  gas  and 
soot  radiative  interchange,  and  absorption  of  soot  radiation  by  the  fuel. 
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Abstract — An  expression  is  derived  for  the  radiative  source  term  governing  the  interaction  of 
molecular  gas  tend  radiation  and  flow  in  nonhomogeneous,  plane-parallel  reacting  flow 
problems.  The  divergence  of  the  net  radiative  flux  is  formulated  in  terms  of  wide-band 
absorptance  model  parameters  for  combustion  products,  and  is  valid  for  all  degrees  of  optical 
thickness.  When  optical  thickness  is  finite,  the  net  absorption  is  obtained  by  integrating  the 
radiation  field  solution  over  the  band  (ineshapes  and  taking  hemispherical  averages.  Illustra¬ 
tive  calculations  for  counterflow  diffusion  flames  will  be  discussed. 


INTRODUCTION 

A  recognized  effect  of  radiation  from  flames  is  that  substantial  fractions  of  energy  can  be  lost  from 
the  flame  environment.  In  addition  to  radiative  loading  of  enclosure  walls,  this  energy  loss  leads 
to  a  temperature  reduction  and  therefore  changes  in  local  gas  density  and  flow  velocity  as  well. 
In  extreme  cases,  the  reduction  in  flame  velocity  can  lead  to  a  substantial  reduction  in  flame  length. 
In  addition,  the  temperature  reduction  affects  flame  chemistry  due  to  the  Arrhenius-exponential 
dependence  on  temperature.  Errors  in  temperature  could  lead  to  substantial  inaccuracies  in 
calculating  formation  rates  of  nitric  oxide,  for  example.  While  all  of  these  effects  can  be  expected 
to  be  most  extreme  in  sooting  flames,  non-luminous  gas  band  radiation  can  be  important  as  well,' 
and  the  way  in  which  one  couples  the  radiation  and  flow  calculations  for  quantitative  assessment 
of  these  effects  in  plane-parallel  or  boundary  layers  of  combustion  gas  products  is  the  subject  of 
this  paper. 


ANALYSIS 

The  coupling  between  radiation  and  flow  occurs  mathematically  in  the  flow  energy  equation,  and 
can  be  represented  in  all  generality  in  the  form 


DH0  _3qf 
D  t  dx, 


(1) 


where  DHJDt  is  a  generalized  representation  of  the  rate  of  change  of  total  enthalpy  (H0)  which 
can  include  convection,  transport,  and  chemical  enthalpy  release,  and  dq?/dxt  is  the  divergence  of 
the  net  radiative  flux.  From  the  equation  of  radiative  transfer  it  is  possible  to  show  that 

=  £  do>  |  dfl*,(a>)[/b<(u)-/(n,<u)] 

/b(a>)=  1.1925  x  !0-5tuJ/(exp(1.438a>/r)-  1)  (c.g.s.;  co  in  cm'1)  (2) 


where  Kt  is  the  local  absorption  coefficient,  7b  is  the  local  Planck  body  function,  and  /  is  the  local 
radiative  flux.  The  integrals  over  all  frequencies  and  directions  of  radiation  propagation  ((2)  must 
be  taken.  In  the  optically  thin  limit,  the  second  term  can  be  neglected,  and  the  isotropic  source 
term  is  equivalent  to 

*4*  j  da>Kt(o))Ib{(o).  (3) 
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Only  frequency-integrated  quantities  are  required,  and  thus  wide-band  models  are  used  for  the 
absorption  bands  of  the  species  that  are  important  in  combustion  product  radiation  (C02,  H20, 
and  CO).  In  the  wide-band  formulation  the  local  absorption  coefficient  for  band  j  belonging  to 
species  i  is  represented  by 


(4) 


where  cc0  is  the  integrated  band  intensity  [units:  cm~'/(g  cm'2)],  Acoy  is  the  band-width;  p,  is  the 
mass  density  of  optically  active  chemical  species  i;  c off  is  the  band  center;  and  C0  =  1  or  2  depending 
on  whether  the  band  is  one-sided  or  symmetric.  With  the  assumption  that  /b  varies  more  slowly 
with  frequency  than  the  absorption  bandshape,  Eq.(3)  can  be  simply  expressed  as:2 


dx, 


=  4«  £  b#» 

EMISS  v 


(5) 


and  the  summation  is  over  all  bands  and  species,  and  is  the  Planck  function  evaluated  at  the 
band  center.  Tabulations  of  the  band  model  parameters  for  the  important  vibrational  infrared  (i.r.) 
and  pure  rotational  transitions  and  the  way  in  which  one  computes  their  temperature  dependences 
are  given  in  Refs.  3  and  4.  We  have  nominally  included  in  our  model  all  15  bands  listed  for  H20, 
C02,  and  CO.  The  optically  thin  or  Planck  limit  expression  given  in  Eq.  (5)  is  convenient  because 
it  depends  only  on  local  parameters,  and  does  not  depend  on  having  an  optical  field  solution.  The 
coupling  between  radiation  and  flow  in  counterflow  diffusion  flames  has  been  examined  in  this 
emission-dominated  limit  in  Refs.  5  and  6,  and  it  has  been  found  that  radiation  can  severely  depress 
peak  temperatures  and  have  a  strong  influence  on  NO,  production  for  low  strain  rates.  However, 
very  low  strain  rates  correspond  to  thicker  flames  where  optical  thickness  effects  could  become 
important,  and  the  purpose  of  the  rest  of  this  paper  is  to  derive  an  expression  for  the  absorption 
term  in  Eq.  (2)  for  plane  parallel  layers  where  all  medium  properties  vary  in  one  coordinate  only, 
say  y: 

=  -  (dco  I  df!A;(cu)/(ft,a>).  (6) 
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The  way  in  which  this  is  done  is  to  derive  a  solution  for  the  frequency-dependent  optical  field,  carry 
out  the  hemispherical  average  by  interchanging  orders  of  integration,  and  then  perform  the 
frequency  integral  in  Eq.  (6).  The  coordinate  system  employed  is  given  in  Fig.  1;  a  boundary  layer 
problem  in  which  medium  properties  vary  only  in  the  y -direction  is  assumed,  and  cold,  black  walls 
will  also  be  assumed  in  the  interests  of  simplicity.  The  effect  of  hot  walls  will  be  discussed  later. 
It  is  convenient  to  carry  out  the  calculation  for  a  single  band,  and  then  sum  over  all  bands/species 
at  the  end;  a  further  approximation  made  here  is  to  neglect  band  overlap  effects  which  are  usually 
relatively  small.  At  a  given  depth  y,  the  sum  of  the  intensities  with  direction  cosines  p  and  —p, 
respectively,  can  be  formally  represented  by 


l„(p.aj)  +  I, p,<a) 


-l 


/bO^to) 


dty  d y' 

d (y'ln)  P 


■r 


-d^  d y' 

d {y'lp)  p 


(7) 


where  r„  is  the  transmissivity  given  for  wide-band  models  by 
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■-if. 
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Fig.  1.  Coordinate  system  for  radiation  analysis  and  opposed  jet  Same  geometry. 


where  fu  and  Pij  are  the  equivalent  line  width  parameter  and  equivalent  broadening  pressure  for 
the  band,  respectively.  The  band  model  parameters  are  here  represented  by  overbars  to  indicate 
that  they  are  equivalent  homogeneous  path  parameters  whose  method  of  calculation  from  mean 
properties  will  be  discussed  later.  Modak7  has  shown  that  it  is  a  good  approximation  in  flame 
calculations  to  assume  the  high  broadening  limit  for  which 

r#  =  exp^-^(<u)^-~ (9) 

This  approximation  should  be  particularly  good  for  the  high-pressure,  counterflow  flame  calcu¬ 
lations  to  be  given  later.  In  this  limit,  Eq.  (7)  is 

/Oi,cu)  +  /(-^,£o)  =  j*  dy7b(y\w)£(co)exp(— K,(a>)|y -y'\lp)dy'lp  (10) 

where  the  ij  subscripts  have  temporarily  been  suppressed  in  the  interests  of  notational  simplicity. 
The  orientation  average  in  Eq.  (6)  is  equivalent  to  the  hemispherical  average 

£  dft{Eq.  (10)]  =  2n  £  df<[Eq.  (10)].  (1 1) 

Performing  this  integration  and  making  the  transformation  z  =  exp(-C0|a>  -  w"”| /Aw),  for  the 
frequency  integration,  Eq.  (6)  is  equal  to 

l  *  » -y'2)  il2) 

where  £,  is  the  exponential  integral  and  the  slowly  varying  Planck  function  has  been  moved  outside 
the  z-integral.  The  z -integral  is  given  by‘ 


[y(l  +  AcD/Acu,  u)lu°  +  £,(«)] 


(13) 
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where  u  =  dp/Ad3[y  -  y'\  is  the  band  center  optical  depth  and  y  is  the  incomplete  gamma 
function  as  given  in  Abramovitz  and  Stegun.9  It  is  more  convenient  to  use  the  y*  form  given 
for  which9 


y(l  +  Adj/Ato,  u)/k(I  +  A‘S/A‘”>  =  f(l  +  Adi/Act))y*(l  +  Adi/Act),  u).  (14) 


The  function  7  *  function  can  be  efficiently  calculated  using  the  series  representation  of  Ref.  9. 
Summing  over  all  species  and  bands  gives  the  final  result  for  the  net  absorption  of  radiant  energy 
by  the  medium 


djl 

8y 


d y'i*A*p)y 


x  [f(l  +  Ad>/Aa>)y*(l  +  Adi/Aeo,  u)  +  £|(u)]#(l  4-  Ad>/Aa))j~‘ 


(15) 


where  u  =  u,y  =  (ap/Aa>)s|y  -y'\. 

The  calculation  of  the  equivalent  homogeneous  path  band  model  parameters  (a,  p,  and  Adi)  is 
based  on  the  Curtis-Godson  approximation  as  given  by  Felslce  and  Tien.10  The  temperature  and 
species  densities  are  simple  path  averages  between  point  y  and  the  point  y'  in  the  running  integral; 
a  =  a(T),  and  Adi  =  £Aa>(T),  where  if  is  a  correction  factor  (Refs.  3  and  7)  ranging  from  1.0  to 
1.44.  The  procedure  for  calculating  if  by  bilinear  interpolation  based  on  the  nominal,  band-center 
optical  depth  and  equivalent  line  width  parameter  (l(T)  is  given  in  Refs.  3  and  7.  For  the  pure 
rotational  band  of  H20  the  band  model  parameters  of  Ref.  1 1  are  used;  for  this  band,  <f  =  1.  The 
calculated  results  to  follow  were  not  particularly  sensitive,  however,  to  the  inclusion  of  if. 

An  alternative  form  of  Eq.  (15)  is  derived  by  converting  the  integration  over  y’  in  Eq.  (7)  to 
an  integration  over  t.  If  one  has  medium  properties  on  a  grid  of  points  denoted  by  the  index  1, 
Eq.  (15)  can  be  shown  to  be  equivalent  to 

=  -ZKlfeoo)  I  Adi  fT§(Yv(U^y)-Yl,(L-*y)) 

ABS  ti  \a(0  /V  I 

-2*1  ^(3'))jAt3‘'>/htf(rwl)r(,(l  -*y)  +  Aw^I^T^r^y))  (16) 

where 


r  =  r(Ad3/Aaj)y*(Adi/Atu,  u)  -  - 

1  4-  Aa)/Ato 


(£,(«)  +  r(l  +  Ad>/Acu)v*(l  +  Aco/Aco,  a)) 


and  it  has  been  assumed  that  all  path  averages  are  taken  from  point  y  to  the  grid  midpoints,  and 
that  Acu  and  /b  are  evaluated  at  the  midpoints.  The  functions  YJ, U+-*y)  and  !*  (/_  — *jf)  are  based 
on  which  of  the  grid  point  pairs  is  closest  to  and  farthest  from  y ,  respectively.  The  effect  of 
absorption  of  hot,  black  wall  radiation  has  been  introduced  in  the  terms  involving  the  wall 
temperatures  Twl  and  Twl;  it  is  understood  that  the  path-averaged  properties  which  go  into  Adi  and 
Y  in  these  terms  are  those  along  the  paths  from  each  wall  to  point  y. 

Equations  (15)  and  (16)  can  be  time-consuming  to  use  because  of  the  running  integrals  that  must 
be  evaluated  for  each  point  y.  The  “analytic”  approximation  of  Modak7  provides  increased 
efficiency  and  reasonable  accuracy  when  the  radiating  species  are  concentrated  in  a  high 
temperature  reaction  sheet  as  encountered  in  diffusion  flames.  Returning  to  Eq.  (7),  it  is  reasoned 
that  most  of  the  contribution  to  the  integrals  will  come  from  regions  near  the  maximum  flame 
temperature  because  of  the  rapid  variation  of  Ib  with  temperature.  The  transmissivity  derivatives 
will  be  more  slowly  varying,  and  can  be  factored  out,  assigned  their  values  at  y  =  y,,  the  point  of 
maximum  flame  temperature.  This  gives,  for  cold,  black  walls: 
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where 
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< 


for  .  ^  if. 
for  t  <  _vf. 


and  the  property  averages  are  calculated  between  y  and  yf.  Equation  (17)  gives  results  nearly 
equivalent  to  Eqs.  (15)  and  (16)  with  much  increased  efficiency,  subject  to  one  restriction  that  will 
be  discussed  in  the  next  section.  Absorption  of  hot  wall  emission  can  be  represented  by  the  same 
terms  as  in  Eq.  (1 1).  These  experiments  are  not  limited  to  slowly  varying  medium  properties  per 
se,  but  would  presumably  have  the  inherent  limitations  of  the  Curtis-Godson  approximation  in 
this  regard. 


SAMPLE  CALCULATED  RESULTS 

Example  calculations  have  been  done  for  an  opposed  jet  diffusion  flame,  as  illustrated  in  Fig.  1 . 
The  configuration  consists  of  two  opposed  fuel  and  oxidizer  jets,  and  is  known  to  lead  to  a 
one-dimensional  structure  normal  to  the  surface  of  the  stoichiometric  mixture.  In  a  small  region 
near  the  stagnation  streamline,  temperature  and  species  concentrations  are  nearly  one-dimensional. 
Temperature  and  radiating  species  concentrations  have  been  computed  using  a  counterflow  flame 
code  developed  by  Smooke.12  The  calculations  are  described  in  more  detail  in  Ref.  13,  where  they 
formed  the  basis  for  a  flamelet  approach  to  predicting  NO,  concentrations  in  an  aircraft  gas  turbine 
combustor.  The  pressure  is  10.5  atm,  and  the  fuel  temperature  has  been  chosen  to  give  an  adiabatic 
flame  temperature  approximating  that  of  jet  fuel.  The  temperature  and  radiating  species  profiles 
are  shown  in  Figs.  2  and  3  for  the  base  case  used  in  these  calculations.  The  latter  corresponded 
to  that  flamelet  in  the  calculations  of  Ref.  13  having  the  lowest  value  of  strain  rate.  The  strain  rate, 
denoted  by  the  symbol  a,  is  the  gradient  of  the  velocity  normal  to  the  flame  structure,  and  is  a 
measure  of  the  rate  at  which  reaction  products  are  pulled  away  in  a  direction  parallel  to  the  reaction 
sheet.  In  the  base  case  shown  in  Figs.  2  and  3,  it  has  the  value  56.1  sec-1.  To  investigate  optical 
thickness  effects,  thicker  flamelets  corresponding  to  lower  strain  rates  were  synthetically  generated 
by  simply  'tretching  the  coordinates  of  the  base  case  by  various  factors.  At  the  low  strain  rates, 
the  species  profiles  correspond  nearly  to  equilibrium  values,  and  it  is  felt  that  this  approximate 
procedure  will  be  adequate  for  illustrative  purposes,  rather  than  using  the  Ref.  12  code  to  get  the 
exact  solution.  Characteristic  dimension  will  scale  with  strain  rate  approximately  as  a-''7.14  No 
coupling  between  radiation  and  the  flame  solution  is  accounted  for  here;  the  energy  equation  in 
the  counterflame  solver  omits  the  radiative  source  term. 
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Fig.  2  Temperature  distribution  for  base  case  Strain  rate  is  56.1  sec*1. 
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Fig.  3.  Radiating  species  profiles  of  base  opposed  jet  case. 

Figure  4  shows  the  total  divergence  of  the  radiative  flux  as  a  function  of  position.  The  position 
is  presented  in  normalized  form  to  correspond  to  the  coordinates  of  the  base  case.  For  the  latter, 
both  the  optically  thin  and  thickness-corrected  solutions  are  given.  Thickness  corrections  are  seen 
to  be  starting  to  be  important  for  this  case.  Two  other  thickness-corrected  solutions  are  given  for 
coordinate  stretchings  of  10  and  25  which  correspond  to  strain  rates  of  0.56  and  0.09  sec-1, 
respectively.  It  can  be  seen  that  thickness  corrections  are  certainly  important  for  strain  rates 
somewhat  lower  than  56.1  sec-1.  It  is  in  this  range  that  significant  temperature  reductions  relative 
to  the  adiabatic  case  are  predicted  with  optically  thin  analysis.5  Optical  thickness  effects  will  cause 
the  temperature  reductions  to  vary  somewhat  less  dramatically  with  decreasing  strain  rate.  Note 
that  there  is  a  small  net  absorption  in  the  relatively  cold  combustion  product  regions  at  the  edges 
of  the  flame  sheet.  In  a  future  publication  the  results  of  incorporating  the  analysis  presented  here 
into  the  counterflow  flame  solver  and  assessing  the  impact  of  radiative  loss  on  NO,  levels  will  be 
addressed,  and  it  will  be  shown  that  very  precise  predictions  require  thickness  corrections  for  low 
strain  rates.15  In  these  calculations  the  water  pure  rotational  band  was  left  out;  it  was  found  that 
with  it  the  factorization  analytic  expression  [Eq.  (17)]  lost  accuracy.  Equation  (5)  is  also  not 
accurate  for  this  band.  Differences  of  a  few  percent  were  found  using  the  exact  solution  [Eq.  (16)] 
with  and  without  this  band. 


-02  -0.1  0.0  0.1  02 
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Fig.  4.  Radiation  source  term  (divergence  of  net  radiative  dux)  as  function  of  strain  rate. 
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Radiative  Dissipation  in  Planar  Gas-Soot  Mixtures 
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Abstract 

The  coupling  of  flame  radiation  and  flow  is  a  subject  of  renewed  interest,  with  an  acceptance 
of  the  fact  that  reacting  flow  models  are  not  complete  without  radiation  effects.  The  interest 
goes  beyond  the  ability  to  predict  heat  transfer  to  enclosure  walls.  It  is  known,  for  example,  that 
substantial  fractions  of  flame  energy  can  be  converted  to  radiation  (Ref.  1),  and  that  the  gas 
cooling  resulting  from  this  non-adiabatic  loss  can  affect  flame  chemistry,  leading  to  inaccuracies  in 
prediction  of  pollutants  like  nitric  oxide.  Neglect  of  radiation  in  flame  hydrodynamic  simulations 
can  also  lead  to  large  errors  in  flow  velocity  and  flame  length.  These  effects  are  expected  to  be  most 
extreme  in  sooting  flames,  but  both  gas  and  soot  radiation  have  to  be  included  in  the  most  general 
model.  Quantitative  assessment  of  these  effects  requires  inclusion  of  an  energy  sink  term  in  the  flow 
energy  equation.  This  sink  term,  the  divergence  of  the  net  radiative  flux,  is  the  difference  of  an 
emission  term  valid  in  the  optically  thin  limit,  and  a  self-absorption  term  that  becomes  important 
when  optical  thickness  effects  need  to  be  accounted  for.  The  expression  for  the  emission  term  is 
generally  trivial  even  when  both  soot  and  gas  band  radiation  are  present,  and  does  not  depend  on 
flame  geometry  because  only  local  properties  are  required.  The  self-absorption  term  is  more  difficult 
because  a  solution  for  the  radiative  flux  is  implicitly  required,  and  integrations  over  all  frequencies 
and  flux  propagation  directions  are  required.  This  term  is  complicated  further  when  both  soot 
and  gas  radiators  are  present  because  of  absorption  profile  overlaps;  the  broadband  soot  absorption 
and  emission  profile  will  overlap  those  of  the  relatively  narrowband  molecular  resonances.  In  this 
paper,  a  complete  solution  for  the  radiative  dissipation  is  given  for  the  one-dimensional,  plane- 
parallel,  or  boundary  layer  flow  problem.  Both  gas  band  and  soot  radiation  are  considered,  and 
the  resulting  solution  is  valid  for  all  degrees  of  optical  thickness.  The  expression  is  given  in  semi- 
analytic  form  in  terms  of  single  quadratures  normal  to  the  flame  structure.  The  theory  will  be 
illustrated  by  example  calculations  for  a  model,  high  pressure,  sooting  counterflow  diffusion  flame. 
Other  applications  might  include  fire  radiation  from  smoke  layers  (Ref.  2),  and  two-dimensional, 
turbulent-jet  diffusion  flames  in  the  boundary  layer  approximation  (Ref.  3). 

Analysis 

The  main  features  of  the  theoretical  analysis  are  as  follows.  The  absorption  terms  are  derived 
by  taking  hemispherical  averages  of  frequency-dependent  solutions  of  the  equation  of  radiative 
transfer,  and  integrating  them  over  the  absorption  bandshapes  of  both  the  molecules  and  the  soot. 
The  soot  absorption  coefficient  is  assumed  to  have  the  Rayleigh  form,  with  scattering  ignored.  For 
nonuniform  paths,  the  soot  terms  can  be  derived  without  approximation;  for  the  gas  band  terms 
a  form  of  the  mean  properties  assumption  is  employed.  The  gas  absorption  profiles  are  assumed 
to  have  the  wideband  model  form,  and  the  transmissivites  are  appropriate  to  the  high  pressure¬ 
broadening  limit.  More  specifically, 

1.  It  is  assumed  that  the  molecular  absorption  bandshapes  are  described  by  wideband  models 
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for  which  (Ref.4)  the  absorption  profile  is  given  by 


K.y  (a/)  = 


QtjPi 

Awjj 


exp 


(1) 


where  ay  is  the  integrated  band  intensity  of  the  jth  resonance  of  radiating  species  i  (HjO,  CO},  CO); 
pi  is  its  mass  density;  Awy  is  the  bandwidth  of  the  resonance,  and  u)(^  its  band  center  frequency,  co 
=  1  or  2  depending  cn  whether  the  band  is  symmetric  or  one-sided.  It  will  further  be  assumed  that 
the  absorption  profiles  vary  much  more  rapidly  with  frequency  than  other  characteristic  functions 
such  as  Planck  function  or  soot  absorption  coefficient,  so  that 

c0  Aw*1  dw  exp  f(u>)  =  f(«40))  (2) 


The  molecular  transmissivity  r  is  taken  to  have  the  form  appropriate  to  the  high  pressure¬ 
broadening  limit  for  which  the  integrated  absorptance  Ay  =  /  dw  (1  -  ry)  has  the  form  (Ref.  5) 


where 
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(3) 


y  being  the  optical  path  length,  Ei  the  exponential  integral,  and  7  the  Euler-Mascheroni  constant. 
Its  derivative  with  respect  to  y  is  thus 
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(4) 


The  suitability  of  this  approximation  for  combustion  radiation  calculations  has  been  discussed 
in  Ref.  5.  It  has  been  found  to  give  good  agreement  in  most  cases  with  more  elaborate  calculations 
based  on  narrowband  models  even  at  one  atmosphere.  Many  of  the  interesting  applications  will 
be  for  high  pressure  systems,  and  the  approximation  should  be  particularly  good  for  these  For 
non-homogeneous  optical  paths,  evaluation  of  the  band  parameters  which  appear  in  Eqs.  3  k 
4,  are  based  on  simple,  path-averaged  temperature  and  density  (Refs.  5,  6).  Overlaps  between 
molecular  bands  are  ignored,  but  the  overlaps  between  the  relatively  narrow  molecular  resonances 
and  the  broadband  soot  absorption-emission  profile  are  included.  It  will  further  be  assumed  that 
the  transmissivity  at  any  frequency  will  be  given  by  the  product  of  the  gas  and  soot  transmissivities, 
e.g.  r(w)  =  r,(w)ry(a/) 

2.  The  soot  absorption  coefficient  i?  taken  to  have  the  form  appropriate  to  Rayleigh  spheroids 


K„s(w)  =  c.wfv  (5) 

where  fT  is  the  particulate  volume  fraction,  the  frequency  in  units  cm-1,  and  the  constant  c,  taken 
to  be  independent  of  frequency.  This  is  equivalent  to  neglecting  the  soot  index  of  refraction.  The 
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Rayleigh  form  of  the  absorption  coefficient  is  suitable  to  first  order  for  aggregate  chains  of  spheroids 
as  well  (Ref.  7).  Light  scattering  by  soot  is  assumed  to  be  small. 

The  generalized  form  of  the  radiative  source  term  is  given  by  the  double  integral  over  wavelength 
and  direction 

f  du  f  dfl  Ks(<t*)  (Ib(w)  —  I(^>w))  =  Q«mi*«  —  Qabs  (6) 

JO  J*  • 

where  Ib  is  the  Planck  function,  I(C,w)  is  the  radiation  intensity,  and  K,(w)  =  Ku(w)  +  Ey  K»y(w). 
Consistent  with  the  above  assumptions,  the  first  term  in  Eq.  6  (emission)  can  be  simply  expressed 
as  (Refs.  8,  9) 


4x  ^  ^  Qiij  Pi  Ibij  4x  Cj  Cl  J 


U4dU  ,  rp5 
- - r  C,fv  1 


(7) 


where  the  factorization  assumption,  Eq.  (2),  has  been  used  for  the  molecular  resonances  ij,  the  soot 
index  dispersion  has  been  ignored  as  mentioned  in  the  soot  part,  and  ci  and  cj  are  constants  in  the 
Planck  function 


Ib(y,w)  =  ciws/(exp(~j)  -  1)  (8) 

Thus 

Ibij=ci(40))3  (exp^f/T)  -  I)'1 

(9) 

=  ci(wy0))3X^  exp(-nc2  wJ0)/T) 

n=  1 

In  the  absorptive  analysis,  one  molecular  resonance  ij  will  be  assumed  to  simplify  the  notation, 
and  a  sum  over  all  ij  performed  at  the  end.  The  radiation  field  solution  is  required  as  a  function  of 
frequency  and  flux  propagation  direction  to  calculate  the  self-absorption  part.  Referring  to  Figure 
1 ,  the  sum  of  the  intensities  with  direction  cosines  p  and  -p  with  respect  to  the  axis  y  normal  to 
the  flame  structure  can  be  formally  represented  as 


/y  ^3  r  rco  _ ^ 

-oo  Ib(y  ’w)d(y7M)  +  Jy  d(7 


-dr_  dy' 

//*)  M 


(10) 


The  absorption  part  will  thus  have  terms  proportional  to  K, 


drij 


aiF/»n 1 


‘■•ij  r*afP^>K»ij  nj  d(yVM)  •  r< 
and  K„  ry  These  represent,  respectively,  gas  absorption  of  gas  emission,  gas  absorption  of 

soot  emission,  soot  absorption  of  gas  emission,  and  soot  absorption  of  soot  emission.  Thus  Q»b»  = 
q  +  Qs*  +  Q.«  +  Q»»  Since  the  8“  and  soot  transmissivities  have  the  form,  respectively, 
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ru (<*>)  =  exp(-K,y(w)|y  -  y'| /y)  =  exp(-a z/fi) 


(11) 


r,(oi)  =  exp(— c,u»lv  (y;y')  |y  -  y'j /m)  =  exp(-b/M) 


(12) 


where 


z  =  zy  =  exp(-Co|w  -  Wy0)  j  /Ata>y) 
b  =  by  =  c.wjj05  fv  (y;y')  |y  -  y'| 

My.y')  =  |y  -  y'"l  \f  My'W' 


all  of  the  terau  contributing  to  Q*b«  will  be  proportional  to 


exp(-(az  +  b  )/n) 


(13) 


The  overbars  in  Eqs.  11,  12  denote  path-averaged  quantities  as  per  Refs.  5,  6.  The  average  over 
all  directions  of  propagation  is  given  by  the  hemispheric  integrals 


/  dfl  |I(w,n)|  =2*  f  dfi  [Eq.10] 

J4*  JO 


(14) 


Thus,  the  exponential  integral  Ei  will  be  a  factor  in  all  terms  from  the  relation 


2jt  /  exp(-(az  -I-  b)/n)  —  =  2*  Ei(az  +  b) 
Jo  M 


(15) 


The  frequency  integrations  are  generally  carried  out  using  /0°°  du>  — »  —■  /J 


Starting  from  the  assumed  forms  for  the  absorption  coefficients  and  transmissivities,  and  mak¬ 
ing  use  of  the  foregoing  assumptions,  the  following  expressions  for  the  terms  contributing  to  the 
absorption  part  are  derived. 

Q««  =  2'  £  (f£)„  /  d>'  u(|^,  b)  (16) 


Q„  =  2*  c.  £  ul0)  J_  <ly'  f.(y')  A«'ii(y')  >bi,  (y')  U  -  1,  »,  b)  (17) 
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Q.g  =  2*  c.  f,(y)  £  wi0>  /_  dy'  ASJ,j(y')  M/)  |y  -  y'l"1  [Ej(b)  -  e,(«  +  b)] 


(18) 


where  a=aij  and  b=bij  are  as  given  in  Eqs.  11,  12,  Ej  is  the  exponential  integral  of  second  order, 
and 

Q*»  =  c,  ej  fv(y)  dy'fv(y')  f°°  do,  a; 6  T  e-ac^'T^  E,(az(w)  +  c.Wfr|y  -  y'|) 
J-°°  Jo 

=  2x  c*  Cl  fv(y)  J  ^  dy'  fv(y')  Te(y')  n~6 

Whf  if_Myl_)kl 

V  Q(y')  '  k  Vp„(yO  +  Q(yO/ 


(19) 


-2jt  cj  Ci  fv  (y)  ^2  (w|j0))5  dy'  fv(y')  ASy  (y')  f)  W) 

ij  ’/-00  n’=l 

•  (E,(a+b)-Ei(b)) 


with 


Pn(y,)  =  nc2/T(y') 

Q(y')  =  c,fv  |y  -  y'| 


Qgg  and  Qgg  both  require  the  integral 


U(p,a,b) 


=  fl  dz 

Jo 


zp  Ei(az  +  b) 


(20) 


Seriei  expansions  for  this  integral  are  given  in  the  Appendix.  The  two  terms  for  Q„,  repre¬ 
senting  the  pure  soot  absorption  plus  a  correction  for  gas  attenuation,  respectively,  are  derived  by 
setting 

E,(az  +  b)  =  E,(b)  +  [E,(az  +  b)  -  E,(b))  (21) 


and  assuming  that 


/ 


dw  (Ej(az  +  b)  -  Ei(b)j  ~  £ 


Awjj  (Ei (a  +  b)  -  Ei(b)) 


(22) 
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where  again  a  and  b  are  as  in  Eqs.  11,  12.  Aside  from  the  factorization  assumption  of  Eq.  2,  Eq.  22 
is  the  only  mathematical  approximation  made.  It  is  possible  to  derive  an  exact  series  representation 
of  this  term  for  b  >  a  which  confirms  the  validity  of  this  approximation. 

Thus  all  terms  contributing  to  field  absorption  are  expressed  in  semi-analytic  form  in  terms 
of  single  integrals  over  the  coordinate  normal  to  the  flame  structure.  These  running  integrals 
are,  however  a  computational  burden.  The  economy  of  the  calculation  is  improved  by  restricting 
attention  to  the  strongest  molecular  lines  (HjO  6.3  and  2.7  micron,  COj  4.3  and  2.7  micron,  CO 
2.35  micron),  which  account  for  most  of  the  radiation.  Limiting  the  terms  in  the  Planck  function 
expansion  Eq.(9)  to  n=2  is  usually  an  adequate  approximation  at  high  temperature  as  well.  This 
type  of  analysis  does  not  lend  itself  to  inclusion  of  non-  black  wall  boundary  conditions.  For 
black  walls,  however,  additional  terms  in  the  expression  for  the  absorption  can  be  calculated.  As 
an  example,  for  a  single  wall  at  temperature  Tw,  the  following  term  is  added  to  the  absorption 
expression  for  the  pure  soot  case 


00  5  _ 

2Jrc1c.fv(y)  4!  £  {(Pnw  +  Qw)-6  -  5QW/P£W  •  [£n(l  +  Pn„/Qw)  -  £  ^(p— - )k]}  (23) 

n=l  k=l  K  Fnw+<'*n 

where 

Pnw  =  ncj/Tw 
Qw  =  c.fy  |y  —  ywJ 

The  net  radiative  flux  to  the  wall  is,  for  pure  soot, 


qtt  =  2w  c«cj 


J  dy'  J  dw  u4  e™"  ^  Ej(c,u;fv  |y'  -  y„|)  fv(y') 

2’c-c‘  fZ  iy'  f’(,,)  „§ 


4.  .  V  1  I 

'  Q(y')  1  ^  k  Vp„(y')  +  Q(>') 


0)' 


(24) 


where 

Pn(y')  =  nc2/T(y') 
Q(y')  =  c,fv|y'  -  y*| 


Corrections  to  Lqs.  23  and  24  for  attenuation  by  gas  band  absorption  would  be  derived  by 
proceeding  as  in  Eqs.  21  and  22,  but  with  Ei  replaced  by  E2 

The  gas  band  contributions  to  net  radiative  wall  flux  may  be  expressed  as 


qR  =  2  k 


dy'  Ibij  (yf)  Awjj  (y')  |yw  -  y'|  1  jE3(b)  -  E3(a  +  b))) 


(25) 
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Sample  Calculated  Results 


Sample  calculations  have  been  performed  for  a  model  opposed  jet  diffusion  flame  as  shown 
in  Figure  1.  The  configuration  consists  of  two  opposed  fuel  and  oxidizer  jets,  leading  to  a  one¬ 
dimensional  structure  in  a  region  near  the  stagnation  streamline  where  application  of  this  one- 
dimension  radiation  analysis  is  appropriate.  Methane  fuel  has  been  assumed,  and  the  flame  prop¬ 
erties  were  calculated  using  a  widely  used  opposed  jet  solver  (Ref.10).  The  flame  was  modified 
by  adding  the  expression  for  the  radiative  source  term  derived  here  to  the  one-dimensional  energy 
equation.  The  clean  flame  calculations  are  described  in  more  detail  in  Ref.  11,  where  they  formed 
the  basis  for  a  flamelet  approach  to  predicting  NOx  concentrations  in  an  aircraft  gas  turbine  com¬ 
bustor.  The  fuel  temperature  was  adjusted  to  give  a  peak  flame  temperature  approximating  that 
of  jet  fuel,  and  the  pressure  of  10.5  atm  is  representative  of  such  a  combustor  at  the  cruise  con¬ 
dition.  In  order  to  assess  the  possible  role  of  soot  in  the  radiative  transfer,  a  sooting  region  has 
been  synthetically  added  to  the  clean  flame  solution  as  shown  in  Figures  1  it  2.  Soot  will  tend  to 
form  on  the  fuel  side  of  the  flame,  <*nd  will  build  up  at  the  stagnation  plane  due  to  thermophoretic 
transport  (Refs.  12,  13).  Calculation  of  the  actual  soot  profile  of  course  requires  a  soot  kinetics 
and  transport  model,  but  this  was  not  available  for  this  study.  Thus,  a  simple  picture  of  the  soot 
profile  is  assumed  in  which  it  has  been  assumed  to  be  uniformly  distributed  between  the  stagnation 
plane  and  a  point  on  the  fuel  side  where  nucleation  might  be  expected  to  be  intense.  Following  Ref. 
14,  we  assume  that  nucleation  will  occur  at  an  elemental  carbon  mixture  fraction  of  .1,  where  the 
adiabatic  temperature  is  about  1675  K.  This  results  in  the  soot  profile  located  shown  in  Figure  2; 
the  effective  radiation  temperature  of  the  soot  is  about  1550K,  nearly  900  K  lower  than  the  peak 
gas  temperature.  The  average  soot  volume  fraction  is  then  treated  parametrically.  The  synthetic 
nature  of  this  calculation  must  be  emphasized;  with  methane  fuel  at  this  temperature  the  flame 
in  reality  would  be  expected  to  be  relatively  clean.  The  intent  here  is  to  gain  some  idea  of  the 
relative  importance  between  soot  and  gas  band  radiation  for  reasonable  ranges  of  the  soot  volume 
fraction.  The  soot  is  coupled  to  the  flame  solution  only  through  temperature  reduction  resulting 
from  the  radiation;  in  reality  the  gas  cooling  due  to  radiation  could  aflfect  the  pyrolysis,  nucleation, 
and  surface  growth  processes  that  grow  the  soot.  An  accurate  treatment  of  this  problem  requires 
that  temperature-dependent  soot  growth  processes  be  added  to  the  chemical  kinetics  mechanism. 
The  basic  gas  kinetics  mechanism  used  here,  where  the  effect  of  radiation  on  NO  concentration  is 
of  interest,  is  a  variation  of  that  proposed  by  Glarborg,  et.  al.,  (Ref.  15).  The  213  reactions  and 
57  species  of  that  mechanism  have  been  culled  to  126  and  30,  respectively.  The  data  set  is  limited 
to  Cl  chemistry,  and  NO2  kinetics  have  been  ignored. 

Opposed  jet  solutions  are  characterized  by  the  strain  rate  (denoted  here  by  the  symbol  a), 
which  is  a  measure  of  the  velocity  gradient  normal  to  the  flame  structure.  It  is  an  indication  of 
how  fast  reaction  products  are  pulled  away  in  the  direction  parallel  to  the  flame  sheet.  High  values 
of  strain  rate  correspond  to  short  gas  residence  times,  thinner  flame  structures,  and  somewhat 
lower  peak  temperatures,  leading  to  reduced  importance  of  radiation.  Low  values  of  strain  rate  are 
characterized  by  the  opposite  attributes,  and  it  will  consequently  be  these  that  are  examined  here. 
In  particular,  a  base  case  value  of  20  sec-1  is  used  for  this  study;  the  adiabatic  temperature  profile 
is  that  shown  in  Figure  2.  For  low  values  of  strain  rate  chemical  equilibrium  radiating  species 
profiles  will  be  a  reasonable  approximation,  and  lower  values  of  strain  rate  can  be  simulated  by 
stretching  the  coordinates  in  proportion  to  a-1,/2.  Figures  3  and  4  are  based  on  fixing  the  adiabatic 
solution  for  a  =  20  sec-1,  stretching  the  coordinates  in  this  way,  and  calculating  the  radiative 
dissipation  according  to  Equations  7  and  16-19.  Figure  3  shows  how  calculated  absorption  and 
emission  vary  with  strain  rate  for  a  clean  flame  condition.  It  can  be  seen  that  optical  thickness 
effects  are  important  even  for  the  a  =  20  sec-1  case,  where  the  characteristic  flame  thickness  is 
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about  1-2  mm.  Figure  4  shows  the  analogous  results  for  a  sooting  flame  where  a  soot  volume 
fraction  of  10*-5*  has  been  assumed.  At  this  value  of  volume  fraction,  the  soot  contributions  start 
to  dominate  the  radiation  process  at  this  pressure. 

Inclusion  of  Equations  7  and  16-19  as  a  power  sink  term  in  the  energy  equation  of  the  opposed 
jet  solver  leads  to  the  predicted  variation  of  radiative  loss  fraction  with  soot  volume  fraction  shown  in 
Figure  5.  As  stated,  soot  radiation  starts  to  dominate  for  volume  fractions  in  excess  of  about  10(~s) , 
with  predicted  loss  fractions  rising  to  about  20%  for  the  highest  values  of  volume  fraction  that  could 
be  reasonably  expected.  Soot  self-  absorption  becomes  important  at  the  highest  volume  fraction. 
Radiative  loss  fraction  is  here  defined  as  the  volume  integrated  net  radiative  power  density  divided 
by  the  volume  integrated  chemical  enthalpy  release  rate.  These  curves  are  of  course  highly  sensitive 
to  the  assumed  thickness  and  effective  radiation  temperature  of  the  soot  zone.  Calculations  at  other 
values  of  strain  rate  are  qualitatively  similar,  with  the  overall  radiative  loss  fraction  decreasing  with 
increasing  strain  rate. 

The  gas  cooling  effects  resulting  from  the  radiative  loss  are  shown  in  Figure  6.  The  percentage 
decrease  in  flame  temperature  will  generally  be  much  less  than  the  radiative  loss  percentage  because 
the  heat  capacity  is  non-linear,  and  because  reductions  in  temperature  will  be  offset  to  some  extent 
by  recombination  of  radicals.  In  the  case  of  soot  radiation,  the  soot  will  be  radiating  from  a 
lower  temperature  region,  and  the  effect  on  the  temperature  there  is  shown  in  Figure  6  for  the 
maximum  value  of  volume  fraction  considered.  Cooling  of  the  flame  zone  will  then  occur  by  thermal 
conduction;  the  percentage  reduction  in  peak  flame  temperature  is  far  less  than  the  20%  radiative 
loss  fraction  for  this  case.  As  stated,  however,  the  temperature  reductions  in  the  soot  zone  are 
probably  large  enough  to  affect  the  soot  growth  processes,  leading  to  a  complicated  coupling  between 
radiation  and  growth. 

The  effect  of  these  temperature  reductions  on  the  profiles  of  NO  is  shown  in  Figure  7.  The 
profile  resulting  from  the  soot  radiative  losses  can  be  contrasted  with  that  resulting  from  assuming 
a  local  non-adiabatic  loss  that  is  everywhere  20%  of  the  local  enthalpy  release  rate.  Because  the 
main  effect  of  the  soot  is  to  cool  regions  where  the  NO  is  not  at  its  peak  levels,  soot  radiation  has 
a  much  smaller  effect  on  NO  reduction  in  these  flames. 
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Appendix 

The  integral 

r1  r1  r°°  e(-**+b)‘ 

U(p,a,b)  =  /  dzzpEi(az  +  b)  =  /  dz  zp  /  dt  -  (Al) 

Jo  Jo  J 1  * 

has  the  following  series  representations 


f  dz  zp  Ei(az  +  b)  =  Ei(a  +  b)/(p  +  1) 

Jo 

(B+b)V2'  /  a  ttn— i(a+  b) 

"Va  +  b/  (p  +  l)n+i 

e~bT(p  +  1)7*(P+  l,a)  _  e~(»+b)  /  b  \  ip/  a  \"  an(a+b) 
(p+1)  (p+ 1)  Vb  +  a/ ^Va  +  b/  (p  +  l)n+i 


(P  >  -1) 


(b  >  a)  (A2) 


(b  <  a) 


where 


(u)n+1  =  T(u  +  n  +  l)/r(u) 
a„(u)  =  n!(l  +  u  +  ua/2!  + . un/n!) 

and  T  and  7*  are  the  gamma  function  and  a  form  of  the  usual  incomplete  gamma  function,  7, 
respectively. 

The  first  series  is  derived  by  performing  the  z-integration  in  Eq.  A1  first,  yielding  (Ref.  16) 


Ai  (p,  at,  1) 


7(p  +  l,a+  1) 
(at)p+1 


(A3) 


7  is  now  expanded  in  a  power  series  (Ref.  17),  and  the  t-integral  is  done  term-by-term.  The  second 
series  is  derived  by  using  Eq.  A3  and  integration  by  parts,  followed  by  term-by-term  integration 
over  t. 
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Figure  1. 
Figure  2. 

Figure  3. 
Figure  4. 
Figure  5. 

Figure  6. 
Figure  7. 


Figure  Captions 

Coordinate  system  for  radiation  analysis  and  opposed  jet  flame  geometry. 

Adiabatic  temperature  profile  in  model  opposed  jet  and  assumed  location  of  sooting  region. 
Strain  rate  =  20  sec-1. 

Variation  of  radiative  dissipation  profiles  with  strain  rate  for  non-sooting  case. 

Variation  of  radiative  dissipation  profiles  with  strain  rate  for  soot  loading  of  lO*-6*. 

Variation  of  predicted  radiative  loss  fraction  with  soot  loading  for  optically  thin  and  thickness- 
corrected  cases. 

Temperature  profiles  for  adiabatic,  gas  band  radiation,  and  gas  band  plus  soot  radiation. 
Predicted  effect  of  radiative  cooling  on  NO  profiles. 
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PREDICTIONS  OF  SOOT  PARTICLE  GROWTH  BASED  ON 
AEROSOL  DYNAMICS  MODELLING 


R.J.  HALL  and  M  B.  COLKET 

United  Technologies  Research  Center 
East  Hartford,  CT  USA 

ABSTRACT 

An  analytical  model  of  soot  formation  in  laminar,  premixed  flames  is  presented  which  is  based 
on  coupling  the  output  of  flame  chemical  kinetics  simulations  with  a  sectional  aerosol  dynamics 
algorithm  for  spheroid  growth.  A  simplified  particle  inception  model  is  employed  in  which  benzene 
plays  an  important  role.  Surface  growth  has  been  based  on  experimental  measurements  and  ab  ini¬ 
tio  calculations  using  various  possible  mechanisms  for  the  surface  chemistry.  Detailed  comparisons 
have  been  made  with  various  flame  data  by  using  experimental  temperature  profiles  and  calculat¬ 
ing  profiles  of  species  concentrations  needed  for  the  inception  rate  and  surface  growth/oxidation 
calculations. 


MODELLING 

Soot  growth  rates  were  calculated  using  a  variety  of  procedures.  Common  to  all  mechanisms  is  the 
assumption  that  soot  mass  growth  occurs  primarily  due  to  acetylene  addition.  As  an  example,  a 
modification  of  the  Frenklach  and  Wang  (1990)  mechanism  was  attempted  in  which  reversibility  in 
the  acetylene  addition  was  included  to  better  explain  high  temperature  growth  data. 

Table  1:  Proposed  Surface  Growth  Mechanism 


Reactions  Considered 

logio(Af) 

Ef 

logio(Ar) 

Er 

1. 

H4-C(s)— •  C(s)4-Hj 

14.40 

12 

11.6 

7.0 

2. 

H4-C(s)  -  C(s) 

14.34 

- 

17.3 

109. 

3. 

C(s)-*  products  4-  CjHj 

14.48 

62 

- 

- 

4. 

CjHj4-C(s)— •  C(s)CHCH 
C(s)CHCH  -  C'(s)  +  H 

12.30 

4 

13.7 

38 

5. 

10.70 

- 

- 

- 

The  resulting  mechanism  includes  possible  acetylene  elimination  from  the  soot  radical  (analogous 
to  phenyl  radical  decomposition)  and  separates  the  acetylene  addition  process  into  a  reversible 
formation  of  the  radical  adduct  and  a  cyclization  reaction.  Assuming  steady-state  conditions  for  all 
intermediate  species,  the  rate  expression  for  soot  mass  growth  is  alculated  to  be 


dm  _  (k,lH[-f  k.a(C,Hal)k«ksxA 

dt  m'  (k _ , ( H 2 J  4  k*(H]  +  k3)k_«k5  4  k«k8|C,li,] 


where  the  Arrhenius  factors  and  activation  energies  for  rate  constants  are  listed  in  Table  1.  In 
Equation  1,  me  is  the  mass  of  a  carbon  atom,  \  is  a  surface  density  of  C,0ot-H  sites  (w  2.3  x  10,s 
cm'1  (Frenklach  and  Wang,  1990),  and  A  is  the  particle  surface  area.  The  rate  of  soot  particle  mass 
growth  will  be  proportional  to  particle  surface  area  because  typical  flame  conditions  correspond  to 


the  free-molecule  regime. 
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Prediction  of  surface  growth  requires  the  local  temperature  and  concentrations  of  species  such 
as  CjHj,  H-atoms,  and  the  oxidants  O2  and  OH.  These,  together  with  the  important  benzene 
concentration,  were  obtained  from  numerical  simulations  of  the  flames  (Harris  and  Weiner,  1983; 
and  Bockhorn,  et  al,  1983)  for  which  soot  data  are  available.  A  comparison  of  the  specific  surface 
growth  rate  predicted  by  Equation  1  with  experimental  data  is  shown  in  Figure  1. 


Height  Above  Burner  (mm) 

Fig.  1.  Comparison  of  specific  surface  growth  rate  predicted  from 
Table  1  mechanism  with  experimental  data 


It  is  assumed  that  small  mass-  or  pre-inception  spheroids  grow  from  benzene  through  surface  growth 
using  the  previously  discussed  particle  growth  models,  with  the  same  gas  phase  rate  coefficients  used 
in  the  derivation  of  the  particle  growth  rates,  G.  In  the  free-molecule  regime,  the  number  density 
distribution  function  n(m,  t)  for  these  spheroids  will  be  given  by  the  classical  growth  law 

^  +  ^(G(t)rn*/3n)  =  ii(m,  t)  (2) 

where  the  r.h.s.  is  the  rate  of  gas  phase  or  chemical  production.  The  solution  to  Equation  2  is 
readily  obtained  using  the  method  of  characteristics.  At  a  mass  threshold  mjnc,  these  pre-inception 
spheroids  are  allowed  to  dimerize/coalesce.  The  inception  rate  is  not  a  well-defined  concept,  but  it 
can  be  considered  to  correspond  to  the  particle  current  at  this  theshold  mass,  e  g. 


on  a  mass  basis.  If  the  chemical  production  term  on  the  r.h.s  of  Equation  2  is  taken  to  be  that 
of  benzene  and  if  it  and  the  surface  growth  rate  do  not  vary  too  rapidly  with  time,  the  solution  to 
Equations  2  and  3  can  be  shown  to  be 

Si(t)  ~  FuJ)Sb'"t(T'}  where  r=f  G(t')dt' 

'  '  G(r  )  J 0  (4) 

r  =  r  -  3(m„,f  -  mb,nt) 

and  nribent  and  Sb«n*  are  the  benzene  molecule  mass  and  mass  source  rate,  respectively.  The 
inception  mass  threshold  is  not  currently  known,  and  is  treated  as  a  parameter  in  these  calculations. 
There  is  considerable  sensitivity  of  the  predictions  to  the  value  chosen  because  it  influences  the  effect 
of  coalescence  on  surface  area,  and  thus  the  rate  of  suspended  mass  growth. 
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The  growth  of  soot  spheroids  has  been  modeled  as  an  aerosol  dynamics  problem,  involving  the  divi¬ 
sion  of  the  size  range  of  interest  (approximately  1  to  100  nm)  into  discrete  intervals  or  classes,  and 
then  solving  a  master  equation  for  the  size  class  mass  densities  with  terms  representing  inception, 
surface  growth  (or  oxidation),  and  coagulation  (coalescence).  Because  the  flames  simulated  corre¬ 
spond  to  free-molecule  surface  growth  and  coalescence,  the  calculation  of  the  sectional  coefficients 
is  much  simplified.  The  dynamical  equations  for  the  sectional  mass  densities  Q t  are  (Gelbard  and 
Seinfeld,  1980;  Gelbard,  1982) 


i=l  j=l  i  =  l 

-  i  Q<  -  Q«  E  *A« Qi 

£  i=£+l 

+  lGtCU-  1  -PtQt 

p  dt 

+  aGi+tN  Qt+tN  -3  Gt  Q<  +  5«Si(t) 


£-1  £-1 


*n  =  -1 
=  +l 


G(t)  >  0 
G(t)  <  0 


where  the  /?  represent  sectional  coagulation  coefficients,  1 G £  are  intra-sectional  growth  coefficients, 
and  *Gfc  are  inter-sectional  growth  coefficients.  Implicit  is  the  assumption  of  complete  particle 
coalescence  after  collision,  with  no  aggregate  or  chain  formation.  With  the  general  functional  form 
for  the  surface  growth,  Equation  1,  the  surface  growth  coefficients  can  be  calculated  analytically  as 


*G£  =  3(fn(mt/mt-i))‘ 


■Gt=-  G.  _ 

rru_tN  —  m£ 


where  the  m£  are  sectional  mass  boundaries,  and  the  m£  represent  the  average  intra-sectional 
masses.  Similarly,  the  coalescence  coefficients  need  to  be  calculated  only  once  to  obtain  the  mass- 
dependencies,  and  then  can  be  scaled  as  \/T(t). 

With  the  solution  for  the  Q«,  a  number  of  aerosol  properties  can  be  evaluated.  The  soot  volume 
fraction,  for  example,  will  be  given  by  fv  =  Q £  where  p,  is  the  soot  specific  density,  assumed 

to  be  that  of  solid  carbon,  1 .8  g/cc.  The  M AEROS  formalism  has  an  extraordinary  ability  to  predict 
volume  fraction  (total  soot  mass)  with  a  relatively  small  number  of  size  classes.  Usually  3-5  are 
enough  to  provide  a  converged  solution,  but  useful  results  can  sometimes  be  obtained  with  as  few 
as  two! 


COMPARISONS  TO  EXPERIMENTAL  DATA 

To  compare  theory  with  the  soot  growth  data  (Harris  and  Weiner,  1983;  and  Bockhorn,  et  al, 
1983),  profiles  of  temperature,  benzene  and  acetylene  concentrations,  and  net  surface  growth  rate 
are  provided  to  the  aerosol  code  as  a  function  of  time  or  height  above  the  burner  surface.  The 
net  surface  growth  rate  consists  of  the  mass  addition  rate  for  acetylene  vapor  deposition,  minus 
oxidative  terms  due  to  oxygen  molecules  and  OH  radicals.  The  latter  typically  are  dominant  low 
in  the  flame. 

For  comparison  with  the  Harris  flames,  CHEMKIN  simulations  (Kee,  et  al,  1985)  were  carried  out 
for  C/O  ratios  of  0.8,  0  92,  and  0.96  since  temperature  profiles  were  available  for  these  flames. 
Experimental  volume  fraction  data  were  available  at  values  of  0.8,  0.84,  0.90,  and  0.94.  The  com¬ 
parison  of  the  Table  1  growth  model  with  the  data  is  shown  in  Figure  2.  It  can  be  seen  that  the 
stoichiometric  dependence  of  the  soot  volume  fraction  is  approximately  satisfied.  The  model  cal¬ 
culations  tend  to  suffer  from  overshoot  at  early  times,  particularly  in  the  C/O=0.8  case,  probably 
because  the  inception  rate  resulting  from  the  simplified  model  is  too  high. 


Fig.  2.  Comparison  of  predicted  soot  volume  fractions  based  on 
Table  1  surface  growth  mechanism  with 
experimental  data  from  Harris  and  Weiner,  1983. 

Comparisons  with  the  high  temperature  acetylene  and  propane  flames  of  Bockhorn,  et  al  (1983)  are 
also  discussed.  It  is  shown  that  agreement  with  the  data  here  is  much  improved  by  inclusion  of  an 
“ageing”  of  the  C,oot  -  H  active  sites  according  to 

~  =  <xp  (-  ['  k.,.(T(t'))dt')  (7) 
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